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Abstract 

o[ 

, Starting from first principle many-body quantum dynamics, we show that the dynam- 

ics of Bose-Einstein condensates can be approximated by the time-dependent nonlinear 
Gross-Pitaevskii equation, giving a bound on the rate of the convergence. Initial data 
are constructed on the bosonic Fock space applying an appropriate Bogoliubov transfor- 
mation on a coherent state with expected number of particles N. The Bogoliubov trans- 
formation plays a crucial role; it produces the correct microscopic correlations among 
the particles. Our analysis shows that, on the level of the one particle reduced density, 
the form of the initial data is preserved by the many-body evolution, up to a small error 
which vanishes as iV -1 / 2 in the limit of large N. 
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A Bose-Einstein condensate is a state of matter of a gas of bosons where a macroscopic 



1 Introduction and main results 

CN 
> 

. fraction of the particles occupy the same one-particle state. The existence of Bose-Einstein 

condensation at small temperature was first predicted in 1925 by Bose and Einstein, who 
considered gases of non-interacting bosons. Seventy years later, in 1995, the existence of 
Bose-Einstein condensates was then confirmed by experiments; see [HE]- Since then, Bose- 
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| Einstein condensates have attracted a lot of attention in theoretical and in experimental 

physics. In particular, they have been used to explore fundamental questions in quantum 
mechanics, such as the emergence of interference, decoherence, superfluidity and quantized 
vortices. In experiments, condensates are initially trapped by strong magnetic fields and 
cooled down at extremely low temperatures (in the nano-kelvin scale). Then, the traps are 
released and the subsequent evolution of the condensate is observed. The goal of this paper is 
to study the dynamics of initially trapped Bose-Einstein condensates. In particular, starting 
from many-body quantum dynamics, we show rigorously that the evolution of the condensate 
can be described, in certain regimes, by the time-dependent Gross-Pitaevskii equation. 

The model. We consider a trapped gas of N bosons, described by the Hamilton operator 

N N 

HT = E (" A ** + V **l?i)) +J2N 2 V(N(xi - x 3 )) (1) 

j=l i<j 
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acting on the Hilbert space L 2 (M. 3N , dx), the subspace of L 2 (M. 3N , dx) consisting of all func- 
tions symmetric with respect to permutations of the N particles. The external potential 
V ex t confines the particles inside the trap. The interaction potential V is assumed to be 
non-negative, spherically symmetric and decaying sufficiently fast at infinity. We denote by 
ao the scattering length of the potential V. To define the scattering length, we consider the 
solution of the zero energy scattering equation 

-A + M/ = (2) 

with the boundary condition f(x) — > 1 as \x\ —> oo. The scattering length is then given by 

8vra = J dxV(x)f(x). (3) 

It is easy to check that, for large \x\, 

f( x ) = l-^ + 0(\x\- 2 ). (4) 
\x\ 

By simple scaling, we find then 

N 2 \ 
-A + — V(Nx) J f(Nx) = 

which implies that the scattering length of the rescaled potential N 2 V(N.) appearing in dTJ 
is given by a = oq/N. 

Ground state properties. Let 

E N = min (ip N , H^iPn) 
Un\\=1 

denote the ground state energy of ([T]). It was proven in [19] that 

lim — — = min £Qp(<f) 
llip||a=l 

where 

£ G p(<f) = j dx [\V<p(x)\ 2 + V ext (x)\<p(x)\ 2 + 4vra |^(x)| 4 ] (5) 

is the Gross-Pitaevskii energy functional. Hence, in the leading order, the ground state 
energy per particle depends on the interaction potential V only through its scattering length 
ao- In [18], it was also shown that the ground state of ([T|) exhibits complete Bose-Einstein 
condensation in the minimizer of the Gross-Pitaevskii energy functional ([S}, in the sense 
that 

7$ -> |^gp)(0gp| (6) 

where |</>gp)(</ ) Gp| is the orthogonal projection onto the (normalized) minimizer 0qp £ 
L 2 (M 3 ) of ([5]) and where 7^ denotes the one-particle reduced density associated with the 
ground state i/jn G L^(R 3Ar ) of (pTJ) , which is defined as the non- negative trace class operator 
with integral kernel 

7^ ) (x;x'):= dx2...dx N ip N (x,x 2 ,...,x N )il; N (x',X2,...,x N ). (7) 
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We assume here that [|^jv|| = lj an d therefore that TV 7^ = 1. The convergence in ([6]), 
which hold in the trace-class topology, implies that, in the ground state of dJJ), all particles, 
up to a fraction which vanishes in the limit of large N, are in the same one particle state, 
described by the orbital <^>gp- The results of [W\ I18| show that the Gross-Pitaevskii theory 
correctly describes the ground state properties of the Hamiltonian (TTJ). 

Time evolution. When the traps are switched off, the system starts to evolve, the dy- 
namics being governed by the iV-particle Schrodinger equation 

id t ipN,t = H N tp Nt t (8) 

with the translation invariant Hamiltonian 

N N 

H N = Y, "Ax, +^N 2 V(N( Xi - Xj )). (9) 

j=l i<j 

It turns out that the time evolution of an initial data exhibiting complete condensation can 
be described, in the limit of large N, by the Gross-Pitaevskii theory. In fact, the following 
result was established in [51 [7J El El Q2] , and, in a slightly different form, in [21]. Consider a 
family ip^ 6 L^(M 3iV ) with bounded energy per particle 

& N ,H N rp N )<CN 

and exhibiting complete condensation in a one-particle state ip 6 // 1 (M 3 ), in the sense that 
the one-particle reduced density 7^' associated with ^iv (defined as in fl7J)) satisfies 

as N — > 00. Then, the solution ^jv,t of the Schrodinger equation (jSJ) still exhibits complete 
Bose-Einstein condensation, in the sense that the reduced one-particle density 7IJI associated 
with tpN,t satisfies 

7S,Wl^)(^l ( 10 ) 
as N — > 00, where <pt is the solution of the time-dependent Gross-Pitaevskii equation 

id t (pt = -A<pt + 8ira \(pt\ 2 <Pt- (11) 

This result establishes the stability of complete Bose-Einstein condensation with respect to 
the time-evolution, and the fact that the condensate wave function evolves according to (fTI]) , 
It justifies therefore the use of the Gross-Pitaevskii equation (jlip to predict and describe the 
outcome of the experiments discussed above. 

Mean field regime. The method used in [6j [7J [91 [10] to prove (fT0|) relies on techniques 
first developed to understand the mean-field limit of many-body quantum dynamics. This 
regime is achieved when considering the time-evolution 

id t iP N>t = Hfy Nit (12) 

generated by the Hamiltonian 

N N 

= - A*, + V( Xl - Xj ) (13) 

j=l i<j 
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in the limit of large N. Also in this case, under appropriate assumptions on the potential 
V, complete condensation is preserved by the time-evolution. Here, the evolution of the 
condensate wave function is governed by the Hartree equation 



id t <p t = -Acp t + (V * \ip t \ 2 ) (ft ■ 



(14) 



The first rigorous derivation of (j!4[) was obtained in [23] for bounded interaction potential, 
i.e. under the assumption ||V||oo < °°- The basic idea in [23] was to study directly the time- 

(k) /irav 

evolution of the family of reduced densities "f Nt , k = 1,2, ... ,N, defined similarly to (|7]). 
From the Schrodinger equation (|12p . it is easy to derive a hierarchy of N coupled equations, 

(k) 

known as the BBGKY hierarchy, for the evolution of the densities Jjy t - As N — > oo, the 
BBGKY hierarchy converges, at least formally, towards an infinite hierarchy of coupled 
equations, which is solved by products of the solution of the Hartree equation (|14|) . The 
problem of proving the convergence towards the Hartree dynamics essentially reduces to 
showing the uniqueness of the solution of the infinite hierarchy. This general scheme, first 
introduced in [23j for bounded interactions, was later extended to potentials with Coulomb 
singularities in [11] (and in [3], for bosons with relativistic dispersions). In [3 El El [TO] , the 
same strategy was then applied to analyze the dynamics generated by the Hamiltonian ([9]) 
and to obtain a rigorous derivation of the Gross-Pitaevskii equation (|lip . 
Writing ([9]) as 



one can interpret the Gross-Pitaevskii regime as a very singular mean-field limit, where 
the interaction converges, as N — > oo, towards a delta- function. From the point of view of 
physics, however, these two regimes are very different. While in the mean field limit particles 
experience a large number of weak collisions, the evolution generated by ([I]) is characterized 
by very rare and strong interactions (two particles only interact when they are very close, at 
distances of order A^ _1 from each others, which is much smaller than the typical interparticle 
distance, of order iV -1 / 3 ). As a consequence, it turns out that correlations among particles, 
which are negligible in the mean field limit, play a crucial role in the Gross-Pitaevskii regime. 
To explain this point, let us consider the evolution of the one-particle reduced density 7^ 
which is governed by the equation 



where 7^ \ is the two-particle reduced density, defined similarly to ([7]), with the normalization 

(2) 

Tr jjy t = 1. If we assume the initial state to exhibit complete condensation and we accept 

that condensation is preserved by the time evolution, we should expect and 7^ to be 

(2) 

approximately factorized. In 7^ v however, we also want to take into account the correlations 
between the two particles. Describing the correlations through the solution of the zero-energy 
scattering equation ([2]), we are led to the ansatz 




(15) 




(16) 



<y [ N > t (x;x') ~ (p t (x)<p t (x'), 

( 2") 

T N ' t (x 1 ,x 2 ;x' 1 ,x' 2 ) ~ f(N(xi - x 2 ))f(N(x' 1 - x 2 )) <p t (x 1 )<p t (x2)'iPt( x i)<Pt{ x 2)- 



(17) 
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Plugging this into (|16p . we obtain a new self-consistent equation for (ft, given by 



id m = -Atpt + {N 2 {N - l)V(N.)f(N.) * |^| 2 ) ip t . (18) 

From ©, we have N 2 (N - l)V(Nx)f(Nx) ~ N 3 V(Nx)f(Nx) ~ 8vra (5(x) as JV ^ oo; 
therefore, in this limit, <pt must be a solution of the Gross-Pitaevskii equation We 
observe here that the presence of the factor f(N.), describing the correlations among the 
particles, is crucial in this argument to understand the emergence of the scattering length in 
(fTT]) . We conclude that any derivation of the Gross-Pitaevskii equation (fTT]l must take into 
account the singular correlation structure developed by the many body evolution. In fact, 
understanding the correlations and adapting the techniques of [23J [TU H] to deal with them 
(for example, when proving a priori regularity of the limiting densities) was one of the main 
challenges faced in [61 [U [10] . As we will discuss shortly, correlations also play a major 
role in the approach presented in this paper. 

The coherent states approach. A drawback of the approach used in [61 El [HJ El [10] is 
the fact that it does not give any control on the rate of the convergence of the many-body 
quantum evolution towards the limiting dynamics (described by the Hartree or the Gross- 
Pitaevskii equation). The problem of controlling the rate of convergence is by no means 
of purely academic interest. In experiments, the number of particles N is large (N ~ 10 3 
in typical samples of Bose-Einstein condensates), but of course finite. For fixed N, the 
statement (|10p is empty. Only an explicit bound on the difference 7^ — \ipt)((pt\ can tell us 
whether the limiting (Hartree or Gross-Pitaevskii) dynamics is a good approximation for the 
many-body evolution. In [22j, a different approach to the study of the many-body quantum 
dynamics in the mean field regime was developed, starting from ideas introduced in a slightly 
different context in [15\ I13|. To explain this approach, we consider the bosonic Fock space 

T = (^L 2 (M 3n ,dx 1 ...dx n ). 

n>0 

The idea here is that on T we can describe states with variable number of particles. The 
normalized Fock-state ip = {^ n '} n >i contains n particles with probability HV^H 2 - For 
/ G L 2 (IR 3 ) we introduce, as usual, creation and annihilation operators a*(f) and a(f) 
which act on T by creating and annihilating a particle with wave function /, respectively 
(precise definitions and basic properties are given in Section [5]). For x G M 3 , we can also 
introduce operator valued distributions a* and a x , creating and annihilating a particle at x. 
In terms of these distributions, we define the mean field Hamiltonian 

= J dxV x a* x V x a x + J dxdy V(x - y)a* x a* y a y a x . 

The operator T~Ln commutes with the number of particles operator M, whose action on if) = 
n >o £ J 7 is given by (A/V) = n^ n \ Moreover, when restricted to the sector with 
exactly N particles, coincides exactly with the mean field Hamiltonian (|13p introduced 
above. The advantage of working on the Fock space is that we have more freedom in the 
choice of the initial state. We will use this freedom, by choosing coherent initial states. The 
coherent state with wave function / G L 2 (M 3 ) is defined as the vector W(f)Q G where 
Q = {1, 0, . . . } is the Fock vacuum, and W(f) = exp(o*(/) — a(/)) is the Weyl operator with 
wave function /. A simple computation shows that 



w(m = e-^ 2 £ = e -"/!ii/ 2 { 1, /, ^, . . . } 
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Coherent states do not have a fixed number of particles. However, the expected number of 
particles is given by H/Hl- To recover the mean field regime discussed above, we consider 
therefore the time evolution of an initial coherent state W(y/~N<p)£l, for 92 G L 2 (IR 3 ) with 
\\(p\\2 = 1- It turns out that coherent states have especially nice algebraic properties (due to 
the fact that they are eigenvectors of all annihilation operators). Making use of these prop- 
erties, one can show that, for large N, the mean field evolution ^N,t = e~ % '^ t W(\^N(p)^l 
of an initial coherent state will again be approximately coherent, of the form W(\fNipt)£~l, 
where (ft is the solution of the Hartree equation (|14p with initial data (ft=o = ¥>• 

Moreover, it is possible to express the difference between the reduced densities associated 
with the evolved coherent state and the orthogonal projection onto <pt in terms of a so called 
fluctuation dynamics, defined as the two-parameter group of unitary transformations 

U m \t- s) = W%y/Nip t )e- m N l( f- s )w{y/Nip s ). (19) 

Essentially, the problem of bounding the rate of the convergence of the many-body evolution 
towards the Hartree dynamics in the mean field limit reduces to the problem of controlling 
the growth of the expectation of the number of particles operator M with respect to IA (t; s), 
uniformly in N. To obtain such estimates, we observe that the fluctuation dynamics satisfies 
a Schrodinger type equation 

id t K mi (t;s) =Cf(t)K mi (t;s) 
with U mi (s; s) = 1 and with the time dependent generator 

£t(t) = [id t W*{VNvt)]W{^vt) + W\^vt)UfW{^Nvt) 

dxV x a* x V x a x + \ dx(V * \tp t \ 2 )(x)a* x a x + / dxdyVix - y)ip t (x)7p t (y) a* x a y 



+ / dxdyV{x-y)((p t {x)<p t (y)a x a* y + <p t (x)(p t (y)a x a y ) 

= J dxdyV(x - y)a* x (ip t (y)a* y +Tp t (y)a y ) + Tjjjy / dxdyV{x - y)a%a*a y a :i 



(20) 



The terms on the third and fourth line do not commute with the number of particles oper- 
ator (because the number of creation operator does not match the number of annihilation 
operators). As a consequence, the fluctuation dynamics U mi (t; s) does not leave the number 
of particles invariant. This is not surprising since U mi (t; s) describes fluctuations around the 
Hartree evolution, which are expected to grow with time. Still, under the assumption that 
the interaction V contains at most Coulomb singularities, bounds of the form 

U m{ (t; 0)WW mf (f; 0)^> < Ce K ^ (Af + 1) V) 

for the growth of the number of particles operator (and, actually, also for all its power) were 
proven in [22] . As a corollary, estimates of the form 



Tr 



<^ r (21) 



on the rate of convergence towards the Hartree dynamics followed. Here T^l denotes the 
one-particle reduced density associated with the evolved state ^N,t = e~ l ^ Nt W(\/~Nip)fl. 
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The one-particle reduced density associated with a Fock space vector is defined by the 
integral kernel 

r £W) = ^^ (22) 

It is simple to check that, for A-particle states, (|22p coincides with the definition ([7]). The 
analysis of the time evolution of initial coherent states is also useful to study the dynamics 
of initial data with a fixed number of particles N. For iV-particle initial data ipx G L 2 (M? N ) 
obtained projecting down a Fock-state of the form W(^/Nf)^ onto the TV-particle sector, 
bounds of the form (|2ip were established in [2], extending the ideas developed in [22], for 
arbitrary $ £ J 7 with finite number of particles and energy (note that this class of iV-particle 
states include factorized wave functions of the form ip® N , which are obtained choosing ip = 
f2). Using techniques similar to those proposed in [2T], bounds for the rate of convergence 
towards the nonlinear Hartree dynamics were also obtained in [T7j, allowing also for potential 
with strong singularities. A different point of view on the mean field limit was given in 
[12], where the convergence towards the Hartree dynamics was interpreted as a Egorov-type 
theorem. 

Bogoliubov transformations and the Gross-Pitaevskii regime. It seems natural to ask 
whether the coherent state approach introduced in [22] can be used to obtain a rigorous 
derivation of the Gross-Pitaevskii equation (|lip . providing at the same time bounds on the 
rate of the convergence of the form (|2ip . A major difficulty in this program is immediately 
clear. From (|20p . we notice that there are two contributions to the generator j£]\i(t), one 
arising from the derivative of the Weyl operator W*(y/N ipt), the other from the derivative of 
e -%{t-s)H.N _ The second contribution, given by W*(^/N<p t )HffW(VNip), can be computed 
recalling that Weyl operators act as shifts on creation and annihilation operators (see (|36p ). It 
turns out that this contribution contains a term, linear in creation and annihilation operators, 
having the form 

J dx [-Aipt(x) + (V * \(p t \ 2 )(x)ip t (x)] a* + h.c. (23) 

This term is large (of the order N 1 / 2 ) and does not commute with the number of particles 
operator. With such a term in the generator, it would be impossible to show uniform (in N) 
bounds for the growth of the number of particles. In the mean field regime, however, (|23p 
is exactly canceled by the contribution proportional to the derivative of W*(y/Ntpt), which 
contain the term 

-\^N J dx (id t ip t (x))a* x - h.c. = -VW J dx [-Aip t {x) + (V * \(p t \ 2 )(x)ip t (x)] a* - h.c. 

where we used the equation (jT3|) . As a result, the generator Cpf(t) on the r.h.s. of ([73]) 
contains only terms which, at least formally, are either order one, or smaller. 

To adapt this approach to the Gross-Pitaevskii regime, we define the Fock-space Hamil- 
tonian 

U N = J dxV x a* x V x a x + dxdyN 2 V(N(x - y))a* x a* y a y a x (24) 

and, following (|19p . we naively introduce the fluctuation dynamics 

U GP (t- s) = W%VN^ N) )e-^- s ^W(VN<p^) (25) 

where, for technical reasons which will be clear later, we choose ip[ to solve the modified 
Gross-Pitaevskii equation 

id t 4 N \x) = -A<pf\x) + (N 3 V(N.)f(N.) * \^ t N) \ 2 )(x)^ t N \x) (26) 
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where / is the solution of the zero-energy scattering equation ([2]) (we use the notation y>t to 
distinguish the solution of (|26p from the solution <pt of the limiting Gross-Pitaevskii equation 
(jlip ). Since the solution of (|26p can be shown to converge towards the solution of (jlip . with 
an error of the order TV -1 , control of the fluctuations around (|26p also implies control of 
the fluctuations around (fTTj) . Let jC^ denote the generator of (|25|) . given by sum of the 
two contributions (id t W*(VNipl N) ))W(VNipl N) ) and W*(\fNipf* ) )U N W{VN<Ji N) ). As in 
the mean field regime, the term W*(\^ip[ N ^)7iNW(\^Nip[ N ' > ) contains a large contribution, 
linear in creation and annihilation operators, given by 



N / dx 



-^f\x) + (N*V(N.) * \4 N ^)(x)4 N \x) 



a* x + h.c. 



The term (idtW* (y N y>\ ))W(y/N(p\ ), on the other hand, contains the linear term 

-VN J dx[(idtip { t N) (x))a* x + h.c] 

= -VN J dx [-Aifl N \x) + {N^iN^fiN^ip^ix^ip^ix)^ a% - h.c. 

In contrast with the mean field regime discussed above, here, because of the factor /, there 
is no complete cancellation between the two large linear terms. Hence, the generator £^ p (t) 
of (I25p contains a large contribution, linear in the creation and annihilation operators, of the 
form 



N J dx(N 3 V(N.)(l-f(N.))*\^ N) \^(x)^l N \x)a x + h.c.) (27) 

and it seems impossible to obtain uniform (in N) bound on the growth of the number of 
particles w.r.t. (|25p . The reason for this failure is that we are trying to control fluctuations 
around the wrong evolution. When we approximate e~ lt ^ N W(\^Nip)tfj by an evolved coher- 
ent state W(y/Ntpl N) )^, we are completely neglecting the correlation structure developed 
by the many-body evolution. As a result, fluctuations around the coherent approximation 
W(y/Nip[ N ^) are too strong to be bounded uniformly in N. Since correlations are, in first 
order, an effect of two-body interactions, we are going to approximate them using a unitary 
operator of the form 

T(k) = exp Q J dxdy (k(x,y)a* x a* y -k(x,y)a x a y ) 



for an appropriate k E L 2 (M 3 x M 3 ), which will be interpreted as the integral kernel of 
a Hilbert-Schmidt operator (again denoted by k) on L 2 (M 3 ). The operator T(k) acts on 
creation and annihilation operators as a Bogoliubov transformation (see Sect. 12.21 for precise 
definitions and basic properties): 

T*(k)a(f)T(k) = a(ch(fc)(/)) + a*(sh(fc)(/)) 
T*(k)a*(f)T(k) = a*(ch(k)(f)) + a*(sh(k)(f)) 

where ch(fc) and sh(k) are the bounded operators on L 2 (M 3 ) defined by the (absolute con- 
vergent) series 

ch ^) = E?iv(^) n and sh ^) = E7^TTv(^ 

n>0 ^ n>0 ^ 
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where products of k and k have to be understood in the sense of operators. Inspired by the 
analysis of El [TU] , where correlations were successfully described by the solution of the 
zero-energy scattering equation dH), we define the (time-dependent) kernel 

h(x, y) = -Nw(N(x - y))ip { t N) (x)ip { t N) (y) (28) 

where (pi is the solution of the modified Gross-Pitaevskii equation (|26|) . and where 

w(x) = 1 — f(x) 

with / being the solution of fl2J). We will consider an initial state of the form W{y/~N (p)T(ko)ip, 
for ip with bounded number of particles and bounded energy, and we will approximate its 
time evolution by W (y/N <pt)T (kt)ijj , leading therefore to the fluctuation dynamics 

U(t;s) = T%k t )W%^<pt)e~ i{t - s ^ HN W(^N Vs )T(k s ). (29) 

In this way, the approximating dynamics takes into account the correct correlation structure, 
and we have a better chance to bound the fluctuations. In fact, we will show in Section H] 
that it is indeed possible to obtain a uniform (in N) control for the growth of the number of 
particles w.r.t. the dynamics (|29p . We will see in Section [6] that the large linear contribution 
(|27p . which will still appear in the generator of (|29p . is compensated by a contribution arising 
from sandwiching the cubic term in W*(y/~N ip[ N ^)MnW(\^N ip^ ) between T{kt) and T*(kt). 
In fact, after conjugating this expression with the Bogoliubov transformation, some of the 
(cubic) terms will not be in normal order (a creation operator will lie on the right of an 
annihilation operator). By the canonical commutation relations, normal ordering produces 
then terms which are linear in the creation and annihilation operators; some of these terms 
cancel exactly the large contribution (|27p . Other important cancellations will emerge between 
the quadratic and the not normally ordered quartic terms; see Section [6] for the details. 
The control of the growth of the number of particles w.r.t. (|29p will imply convergence of 
the reduced densities associated with the fully evolved Fock state e~ MNt W(VNift)T(kt)ip 
towards the orthogonal projection onto the solution of the Gross-Pitaevskii equation (jlip. 
with a bound on the rate of the convergence. Note that the idea of adding an operator with 
a quadratic expression in creations and annihilation operators in the exponent to the Weyl 
operator generating the coherent state was used in the mean-field regime in [14J to obtain 
norm-approximations of the many-body dynamics. 



The main theorem. We are now ready to state our main result. 



Theorem 1.1. Let tp € -£f 4 (R 3 ) ; with \\<p\\2 = 1- LetH^ be the Hamilton operator defined in 
P^P, with a non-negative and spherically symmetric potential V £ L 1 n L 3 (M 3 , (1 + |x| 6 )dx). 
Let ip £ J- (possibly depending on N ) be such that 



for a constant D > 0. Let denote the one-particle reduced density associated with 

the evolved state e~ l ^ Nt W{\fN(p)T{kQ)%l), as defined in \2S\) . Then there exist constants 
C,ci,C2 > 0, depending only on V,\\<p\\jj& and on the constant D appearing in i^(J\). such 
that 



(30) 



Tr 



N,t 



< 



Cexp(ci exp(c 2 |t|)) 



(31) 
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for allt £l and N E N. Here (ft denotes the solution of the time- dependent Gross-Pitaevskii 
equation 

id t <ft = -Apt + 8ira \(pt\ 2 (pt (32) 
with the initial condition (ft=o = <f- 
Remarks. 

(i) Let us point out that we insert the correct correlation structure in the initial data. 
Our result implies the approximate stability of states of the form W(y/~N ip)T(ko)ip with 
respect to the many-body evolution (in the sense that the evolution of W (y/~N ip)T(ko)ip 
has approximately the same form, just with evolved (ft, up to a small error). It does 
not imply, on the other hand, that the correlation structure is produced by the time- 
evolution. This is in contrast with the results of El El QZ3] , which can also be applied 
to completely factorized initial data. It remains unclear, however, if it is possible 
to obtain convergence with a A -1 / 2 rate (or with any rate) for initial data with no 
correlations (the problem of the creation of correlation was studied in [5]). 

(ii) The time dependence on the r.h.s. of ([31 p deteriorates fast for large t. This however 
is just a consequence of the fact that, in general, high Sobolev norms of the solution 
of (|32p can grow exponentially fast. Assuming a uniform bound for ||y?t||#4, the time 
dependence on the r.h.s. of ([3T]) can be replaced by C exp(K\t\). 

(iii) To simplify a little bit the computations, we did not include an external potential 
in the Hamiltonian ()24p generating the evolution on the Fock space. In contrast to 
[HJ El [10] , the approach presented in this paper can be extended with no additional 
complication to Hamilton operators with external potential. This remark is important 
to describe experiments where the evolution of the condensate is observed after tuning 
the magnetic traps, rather than switching them off. 

(iv) The convergence (|3ip and the fact that the limit is a rank-one projection immedi- 

(k) 

ately implies convergence of the higher order reduced density T N t , associated with the 
evolved state <S N>t = e- inNt W(VNtp)T(k )ip. Similarly to J22J, is defined as the 
non-negative trace class operators on L 2 (M 3fc ) with integral kernel 

r (fc) , v _ (^N,t,a* yi ■ ■ ■ a* yk ax k ■ ■ ■ a xl ^N,t) 

N > t{XU ' • • ' Xk] Vl ' ' ' ' ' Vk) ~ (^ N ,t,M(M - 1) . . . (AT - k + l)*jv )t ) ' 

Following the arguments outlined in Section 2 of [T7], (|3ip implies that, for every k E N, 



Tr 



k 1 ' 2 



Cl 



To obtain bounds for the convergence of the A;-particle reduced density with the same 
N~ 1 / 2 rate as in (|3ip . following the same approach used below to study would 
require to control the growth of higher powers of the number of particle operator with 
respect to the fluctuation dynamics ([29]) . This may be doable, but the analysis becomes 
more involved. 

(v) Theorem 11.11 and the method used in its proof can also be applied to deduce the con- 
vergence towards the Gross-Pitaevskii dynamics for certain initial data with a fixed 
number of particles. In Appendix [Cl we consider initial A-particle states of the form 
P N W(VNip)T(k )i(;, for ip E T satisfying ([3D]), assuming \\P N W(VN<p)T(k )ip\\ » 
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N~ l l 2 for large N (it is explained in Appendix [C] why this is a reasonable con- 
dition). Here Pn denotes the orthogonal projection onto the TV-particle sector of 
T . It remains to be understood which class of iV-particle states can be written as 
P N W{y/Nip)T(k )ip, for a V G ? satisfying 

2 Operators on the Fock space 

The bosonic Fock space over L 2 (R 3 ) is the Hilbert space 

jr = L 2 (R 3 )®° n = C L^(R 3n ), 

n>0 n>l 

with the convention that L 2 (M 3 )® s0 = C. Here L 2 (R 3n ) is the subspace of L 2 (R 3n ) consisting 
of all functions that are symmetric with respect to arbitrary permutations of the n variables. 
Vectors in T are sequences n >o of n-particle wave functions V {n) G L 2 (R 3n ). The 

inner product on T is defined as 

n>0 

= ipiipj + / " " •d^nV'l (^1) • • • >a;n)^2 (^li • • • ,x n ). 

n>l ^ 

On J 7 , we can describe states where the number of particles is not fixed. The Fock space 
vector ip = V* j • • • } describes a coherent superposition of states with different number 

of particles; the n-particle component is described by ifj^ n \ for any n 6 N (the probability 
that the state has exactly n particles is given by \\ 2 ). A state with exactly ./V particles 
is described on J 7 by a sequence n >o where ip^ 1 = for all n ^ N and ip( N > = ^at. 

The vector {1, 0, 0, ... } € J- is called the vacuum and is denoted by £1. States with a fixed 
number of particles are eigenvectors of the number of particles operator, defined by 

(Nip) (n) = mP (n) . 

To define an evolution on J 7 , we introduce the Hamilton operator 

(?W n) =nP^ 

with the n-th sector operator 

n n 

«l B) = E ~ A *i + E N 2 V(N( Xl - x 3 )). 

j=l i<j 

Note that the subscript N in the notation T~Ln is not related with the number of particles 
(since this is not fixed on the Fock space), but only reflects the scaling in the interaction 
potential (of course, at the end, N will also be related with the number of particles in the 
initial Fock state; otherwise, there would be no relation with the regime discussed in Section 
[TJ. Observe that, by definition, the Hamiltonian T~Ln commutes with jV. As a consequence, 
the evolution generated by Hn leaves each n-particle sector invariant. In particular, 

e-^'{0,...,0,^,0,...} = {0,...,0,e-^ViV,0,...} 

where Hn is the iV-particle Hamiltonian defined in ([9]). In this sense, the iV-body dynamics 
is embedded in the Fock space representation. 
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It is very useful to introduce creation and annihilation operators on F. For / G L 2 
the creation operator a*(f) and the annihilation operator a(f) on T are defined as 

1 n 

(a* (f)ip) (n) (xi, . . . , x n ) = —= ^2 fixjW^ (xi, . . . , Xj-i,x j+ i, ...,x n ), 

VU 3=1 

{a{f)^ n \x l ,...,x n ) = V^+l [ dxj{x)^ n+l \x,x 1 ,...,x n ). 



The operators a*(f) and a(f) are unbounded, densely defined and closed. Note that a*(f) 
is linear in /, while a{f) is anti-linear. The creation operator a*(f) is the adjoint of the 
annihilation operator a(f), and they satisfy the canonical commutation relations 

[a(f),a*(g)} = (f,g) L 2 and [a(f),a(g)] = [a*(f), a*(g)] = (33) 

for f,g G L 2 (M 3 ). We also introduce the self-adjoint operator 

0(/) = a*(/) + a(/). 

Although creation and annihilation operators are not bounded, they are bounded with 
respect to the number of particles operator j\f (actually, to its square root). This is the 
content of the next standard lemma (see [22\ for a proof of this well-known result). 

Lemma 2.1. Let f £ L 2 (M 3 ). Then, for any ip G T , 

\\aWW < H/lkllA/- 1 /^!, 

IKWII<ll/l| 2 ||(A/'+l) 1 /Vl|, (34) 
||0WII<2||/|| 2 ||(AA+1) 1 ^||. 

We will make use of operator- valued distributions a* and a x , with x G M 3 , defined so 
that 

a*(f) = J dxf{x)a* x and a(f) = J dx f(x)a x 
for / G L 2 (]R 3 ). The canonical commutation relations assume the form 



[a x , a* y ] = 5(x - y) and [a x , a y ] = [a*, a*] = 0, 

where 5 is the Dirac delta distribution. 

In terms of the operator valued distributions a* and a x , the number of particles operator 
j\f and the Hamilton operator Hn can be written as 



j\f = I dx a* x a x 



and 



Un = j dxV x a* x V x a x + ^ J dxdyN 2 V(N(x - y))a x a* y a y a x . (35) 

The first term in the Hamiltonian is the kinetic energy; since it will play an important role 
in our analysis, we introduce the notation 



K, = J dxV x a*Va x . 
Note that, like Hn, leaves every n-particle sector invariant, and 

n 

(/CV0 W = 5>A^) 

3=1 

for any ip G T . 



12 



2.1 Weyl operators and coherent states 

For / G L 2 (]R 3 ), we define the Weyl operator 

W(f) = e a * (/) - a(/) 

acting on the Fock space T . In the following lemma we collect some well-known properties 
of Weyl operators. 

Lemma 2.2. Let f,g<= L 2 (R 3 ). 

(i) Weyl operators satisfy the Weyl relations 

W(f)W(g) = W(g)W(f)e- 2ilm{f > 9) L 2 = W(f + g)e~ ilm{f ' 9) ^ . 

(ii) The operator W( f) is unitary on T and 

W(f)* = W(f)- 1 = W(-f). 

(iii) We have 

W(fTa(g)W(f) = a(g) + (gJ) L2 and W(f)*a*(g)W(f) = a*(g) + (f,g) L 2. 

(36) 

Using Weyl operators, we construct coherent states on T. For / G L 2 (M 3 ), the coherent 
state with wave function / is defined as W{f)Sl, where = {1,0, . . . } is the vacuum vector 
in J 7 , describing a state with no particles. Since W(f) is unitary, coherent states are always 
normalized. From the canonical commutation relations (|33|) . it follows that 

w(m = e -n/ni/ 2 e^(/)Q = e-"/iii/ 2 y 

n>0 

= e -ll/ll!/2 J 1 f J_ J_ \ 

In particular, coherent states do not have a fixed number of particles. Instead, they are 
given by linear combinations of states with all possible number of particles. From (|36p . the 
expected number of particles in the state W(f)Q is given by 

(W(m,MW(m) = j dx (n,(a* x +J(x))(a x +f(x))n) = \\f\\l 

More precisely, one can show that the number of particles in the coherent state W(f)tt is a 
Poisson random variable with average and variance given by WfW?,. 

Coherent states have particularly nice algebraic properties, which also simplify the study 
of their time-evolution. These properties are a consequence of the fact that coherent states 
are eigenvectors of all annihilation operators; in fact, from (|36p . we find 

a(g)W(m = W(f)(a(g) + (g, /»fi = (g, f)W(f)Q 

for all f,ge L 2 (R 3 ). 
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2.2 Bogoliubov transformations 

For f,g G L 2 (M 3 ), we introduce the notation 

A(f,g) = a(f) + a*(g). (37) 

We observe that 

(A(f,g)y=A(gJ)=A((° J J Q )(f,g)) (38) 



where J : L 2 (IR 3 ) — > L 2 (R 3 ) is the antilinear operator defined by Jf = /. From the canonical 
commutation relations (|33p . we find that the operators A(f,g) satisfy the commutation 
relations 

[A(f 1 ,g 1 ),A*(f 2 ,g 2 )] = ((f 1 ,g 1 ),S(f 2 ,g 2 )) L2(BL2 with S =( J °_ 1 V (39) 



A bounded linear map 6 : L 2 (R S ) © L 2 (M 3 ) -)■ L 2 (M 3 ) L 2 (M 3 ) satisfying 

e(°, fW° f)e (40) 



x J o y v J 

and 

5 = 6*59 (41) 

is called a Bogoliubov transformation. Bogoliubov transformations are linear maps of the 
pairs (f,g) G L 2 (R 3 ) © L 2 (M 3 ) with the property that the operators b*(f) and b(f), defined 
similarly to (|37p by the equation 

6(/) + fe*(5) = A(G(/,p)) 

for any f,g£ L 2 (R 3 ), are still creation and annihilation operators satisfying canonical com- 
mutation relations and being adjoint to each other. It is simple to check that a general 
Bogoliubov transformation can be written as 



for bounded linear maps U, V : L 2 (R 3 ) L 2 (M 3 ) satisfying U*U-V*V = 1 and U*V = V*U 



Here we use the notation U = JUJ for any bounded operator U on L (R ) (the integral 
kernel of U = JU J is given by U(x,y)). 

For a kernel k G L 2 (R 3 x R 3 ) with k(x,y) = k(y,x), we define now the operator 



T{k) = exp Q J dxdy (k(x,y)a* x a* - k{x,y)a x a y ) S j (43) 
acting on the Fock space T . 

Lemma 2.3. Let k G L 2 (R 3 x R 3 ) be symmetric, in the sense that k(x,y) = k(y,x). 
(i) The operator T{k) is unitary on J- and 

T(k)* = T(k)- 1 = T(-k). 
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(ii) For every f,g G L 2 (]R 3 ) ; we have 

T(k)*A(f,g)T(k) = A(@ k (f,g)) (44) 
->• L 2 (IR 3 ) © L 2 (IR 3 ) is the Bogoliubov transformation 



where 9 fc : L 2 (R 3 ) © L 2 
defined by the matrix 



where ch(k),sh(k) : L L 



( ch(k) sh(k) 
k ~ V sh(k) ch(k) 

- L 2 (IR 3 ) are the bounded operators defined by 



sh(k) = £ 



1 



n>0 



(2n + 1)! 



(kk) n k 



where products of k and k have to be understood in the sense of operators. 
(Hi) We decompose 

ch{k) = 1 +p(k), sh(k) = k + r(k), (45) 

where 1 denotes the identity operator on L 2 (M 3 ). Then p(k) and r(k) (and therefore 
sh(k) ) are Hilbert- Schmidt operators, with 



l|p(*)l|2<e« 



|r(fc)|| 2 <e^ 2 , \\sh(k)\\ 2 < 



(46) 



(Here \\p(k)\\2 denotes the L 2 (R 3 xM 3 ) norm of the kernel p(k)(x,y), which agrees with 
the Hilbert- Schmidt norm of the operator p(k)). 



(iv) Suppose now that k £ L' 



is s.t. V\k G L 2 (M 3 x M 3 ). Then, by symme- 



try, also V 2 A; £ L 2 (R 3 x M 3 ) (we use here the notation (W\k){x,y) = X7 x k(x,y) and 
(V 2 fc)(x, y) = V y k(x, y); note that V\k and V 2 k are the integral kernels of the operator 
products V/c and —kV). Moreover 

||Vip(fc)|| 2 ,||Vir(fc)|| 2 <ell fe ll 2 ||Vi(^)|| 2 , 
||V2p(fc)||2,||V 2 r(fc)||2 <e l|fc||2 ||V 2 (fcA;)|| 2 . 

(v) If the kernel k depends on a parameter t (later, it will depend on time), and if deriva- 
tives w.r.t. t are denoted by a dot, we have 



||p«ll2,||r(A0|| 2 <||fe|| 2 e" 



and 



||Vip(fc)|| 2 , ||Vir(AO|| 2 < Ce^ (\\k\\ 2 \\ Vi(fcfc)|| 2 + ||Vi(fcfc)|| 2 + ||Vi(fefe)|| 2 ) , 
||V 2 p(fc)||2, ||V 2 r(fc)|| 2 < Ce^ (||fc||2||V 2 (fcfc)||2 + ||V 2 (fcfc)||2 + l|V 2 (fcI)|| 2 ) . 

Proof, (i) is clear. To prove (ii), we observe that, letting 

B = i j dxdy (k(x,y)a* x a* y -k(x,y)a x a y ) , 
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we have, for any f,g € L" 



e- B A(f,g)e B = A(f,g) + £ dAi^-e"** A(f, g)e x ^ 
= A(f,g)- f dX 1 e- x ' B [B,A(f,g)]e x ' 

JO 



Iterating, we find 
e~ B A(f, g)e B = A(f, g) + £ ad J B <?)) 



+ (-1 



,n+l 



An 



ciA n+1 e- A "+ lB ad^ +1 (yl(/, 5 )) e A "+ lB 

(47) 

where ad^(C) = C] and ad^ +1 (C) = [B,&d B ((J)]. A simple computation shows that 
a d B (A(f,g)) = [B,A(f,g)}- 

and therefore that 

a d j B (A(f,g)) = (-iyA 



-A 



k 
k 



k 
k 



We have 



k 
k 



2 m 



(kk) m 
(kk) r 



and 



k \ 2m+1 _ ( 0_ _ (£±) m £; 
k ) ~ \ (kk) m k 

for every m E N. Inserting all this in (|47p . we obtain (|44p . if we can show that the error 
converges to zero. We claim, more precisely, that the error term on the r.h.s. of (|47|) vanishes, 
as n — > oo, when applied on the domain D{M 1 ^ 2 ). To prove this claim, we start by observing 
that 

|| + \)^e~ XB (M + l)- 1/2 \\ < e |A|l|fe|12 
for every A£l. Assuming for example, that n is odd, (|48p implies that 

/■Al /-An 

dAi / dX 2 ... dA n+1 e- A "+ lB ad" +1 (^(/,9))e A " +lB (AA+l)- 1 / 2 



< 



(48) 



— \\A ((^)(™ +1 )/ 2 /, (^) (n+1)/2 5 ) (A/" + l)- 1 / 2 ! 



(n + 1) 

<(II/I| 2 + Nl2)ell*i^ 



(n+l)/2| 



< c- 



(n + 1)! - (n + 1)! 

which vanishes as n — >■ oo. The case n even can be treated similarly. To prove (|48p . we 
observe that 
d 



dX 



dxdy k(x,y){e XE ! ip , a* a* e XB ip) — I dxdyk(x,y)(e XB ^,a x a y e XB ip) 



-AB 



,-AB„ 



-AB„ 



< 2 dx ||o T e" 



AB„ 



| a* (fc(x, .))e' 



-AB„ 



<2||fc|| 2 ||(AA + l) 1/2 e 



1/2 -AB, ,,2 
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Gronwall's Lemma implies (08]). 
To prove (iii), we notice that 



{kk) n 



< 



v mrh < 

4^ (2n)! 

n>l v ' 



|2n 



n> 



i (2n) 



< e 



INI 



where we used that, by Cauchy-Schwarz, for any two kernels K i, K 2 G 

2 



|^i^ 2 | 



dxdy 



dzK 1 (x,z)K 2 (z,y) 



(49) 



dxdydzi dz 2 K x (x, z\)K\ (x, z 2 )K 2 {z\ , y)-K" 2 (22 , y) 
< J dxdydz l dz 2 \K 1 (x,z l )\ 2 \K 2 (z 2 ,y)\ 2 

= \\ K 2\\ 2 - 

The bounds for r(k) and sh(k) can be proven similarly. 

To show (iv) we write, using the fact that the series for p(k), r{k) and sh(fc) are absolutely 
convergent, 



V lP (k) = Vi(fefc) 
Vir(Jfe) = Vi(fcfc) 



1 

^ (2n P ' 



.n 

00 



^(2n + l)P ' 



Applying (|49|) . we find the desired bounds. The bounds for the derivative V2 can be obtained 
similarly. 

Finally, to show (v), we remark that 



The bound for r(k) can be proven analogously. From the product rule, we also find that 



71=2 



n=l 



< ell fc H 2 (\\k\\ 2 ||Vi(fcfc)|| 2 + ||Vi(fcfc)||2 + ||Vi(k)|| 2 
The other bounds are shown similarly. 



□ 



3 Construction of the fluctuation dynamics 

In this section, we will construct an approximation for the full many-body evolution of an 
initial data of the form W(y/Ntp)T(ko)ip, as considered in Theorem ll.il Our approximation 
will consist of two parts. First of all, the evolution of an approximately coherent state will 
be approximated by a coherent state with an evolved one-particle wave function (later, we 
will take care of the correlation structure). 
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For a given ip G i? 1 (IR 3 ), we define ip\. as the solution of the modified time-dependent 
Gross-Pitaevskii equation 

= -A^) + (iV 3 /(iV.)F(iV.) * |^| 2 ) if™ (50) 

where / denotes the solution of the zero-energy scattering equation ([2]). For technical reason, 
which will become clear later on, it is more convenient for us to work with the solution of the 
modified Gross-Pitaevskii equation (|5U|) , rather than directly with the solution of (|32p . Since 
N 3 f(Nx)V(Nx) -)• 8ira 5(x), the solution (p[ NS> converges towards the solution of (|32p . as 
— t- oo. This is proven, together with other important properties of the solutions of (|32p 
and (f50j) . in the next proposition. 

Proposition 3.1. Let V G L 1 n L 3 (R 3 , (1 + \x\ 6 )dx) be non-negative and spherically sym- 
metric. Let f denote the solution of the zero-energy scattering equation with boundary 
condition f(x) — > 1 as \x\ — > oo. Then, by Lemma \ 3.$\ below, < / < 1 and therefore 
Vf > with Vf G L 1 n L 3 (R 3 , (1 + \x\ 6 )dx). Let ip G H l {M?), with \\ip\\ 2 = 1. 

(%) Well-posedness. There exist unique global solutions p.,^^ G C(R; i^QR 3 )) o/ i/te 
Gross-Pitaevskii equation 132\) and, respectively, of the modified Gross-Pitaevskii equa- 
tion Ii50\) . with initial data <p. These solutions are s.t. ||<£>*||2 = ||y>[ H2 = 1 / or oil 
t G R. Moreover, there exists a constant C > uii£/i 

(JV) 11 ^ ^ 

/or a// t G R. 

(mJ Propagation of higher regularity. If we make the additional assumption that ip G 
H n (M. 3 ), for some integer n > 2, then ipt,<p[ N ^ G H n (M. 3 ) for every t £ I. More- 
over there exist constants C > depending on \\p\\H n and on n, and K > 0, depending 
only on \\<p\\jji and n, with 

\\n\\Hn,\\4 N) \\H»<Ce K W (51) 

for all t G R. 

(Hi) Regularity of time derivatives. Suppose p G -£f 4 (R 3 ). Then there exists a constant 
C > 0, depending on \\p\\jf4, and K > 0, depending only on \\p\\u\, such that 

\\^\\H^h<Ce^ 

for all t G R. 

(iv) Comparison of dynamics. Suppose now p G i/ 2 (R 3 ). Then there exist constants 
C,ci,C2 > 0, depending on \\p\\h 2 ( c 2 actually depends only on \\<p\\jfi ) such that 

n (AO n / Cexp(ci exp(c 2 |i|)) 

\\ft -<Pt\\2< -Tr 

for all t G R. 

The proof of Proposition 13. II can be found in Appendix IA1 

Using the solution (p^ of (f50j) . we are going to approximate the coherent part of the 
evolution. As explained in the introduction, however, this approximation is not good enough. 
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The many-body evolution develops a singular correlation structure, which is completely 
absent in the evolved coherent state. As a consequence, fluctuations around the coherent 
approximation are too strong to be controlled. To solve this problem, we have to produce 
a better approximation of the many-body evolution, in particular an approximation which 
takes into account the short-scale correlation structure. To reach this goal, we are going to 
multiply the Weyl operator W (y/~N ip[ N ^ ) , which generates the coherent approximation to 
the many-body dynamics, by another unitary operator T(k), having the form (j43f) . obtained 
by taking the exponential of a quadratic expression in creation and annihilation operators. 
The kernel k G L 2 (R 3 x R 3 ) has to be chosen so that T(k) creates the correct correlations 
among the particles. Since correlations are, in good approximation, two-body effects, we can 
describe them through the solution / of the zero-energy scattering equation ([2]). We write 

f(x) = 1 - w(x) (52) 

with limui^jx, w(x) = 0. The scattering length of V is defined as 

8-7rao = / dxV(x)f(x). 



Equivalently, ao is given by 

ao = lim u;(x)|x|. 

\x\— >oo 

Note that, if V has compact support inside {x G R 3 : \x\ < R}, then oq < R and w(x) = 
oq/\x\ for \x\ > R. In general, under our assumptions on V, one can prove the following 
properties of the function w. 

Lemma 3.2. Let V G L 1 nL 3 (R 3 , (l + |x| 6 )dx) be spherically symmetric, with V > 0. Denote 
by f the solution of the zero-energy scattering equation (OJ) and let w = 1 — /. Then 

< w{x) < 1 for all x G M 3 . 

Moreover, there is a constant C > such that 

C C 

w(x) < i — : and \Vw(x)\ < -. — ^ . (53) 

w ~ \x\ + 1 1 Wl ~ \x\ 2 + 1 y J 

Proof. Standard arguments show that < f(x) < 1 holds for every x£l 3 (f( x ) ^ 1 follows 
from V > 0, because of the monotonic dependence of / on the potential; see [201 Appendix 
C]). This implies that < w{x) < 1 for all x G M 3 . From the zero energy scattering equation, 
we have —Aw = Vf/2. This implies that 

w(x) = C [dy —±—V{y)f{y) and Vw(x) = C [ dy^^ V{y)f(y) 
J \x-y\ J 

for an appropriate constant C G R. Using \x\ < \x — y\ + \y\, the fact that / < 1, and the 
Hardy-Littlewood-Sobolev inequality, we find 

|(1 + \x\)w(x)\ <C f dy (-r^—r + 1 + r^-r) V(y)f(y) 

J \\ x -y\ \x-y\J 
<c(||v||3/2 + II^IUi + ll|y|^(i/)|| 3 /2) 

and, analogously, 

|(i + \x\ 2 )w(x)\ <c [dy (-r-^ + 1 + t^-u) v (y)f(y) 

J \\x-y\ z \x-y\ z J 
<C(||y||3 + ||^||i + |||y|V|| 3 ). 
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The right hand side of the last two equations is bounded under the assumption V G L 1 n 
L 3 ((l + \x\ 6 )dx). □ 

The zero-energy scattering equation for the rescaled potential N 2 V(Nx) is then solved 
by f(Nx). We define w{Nx) = 1 - f(Nx). Clearly 

lim w(Nx)\x\ = —j, 
|as|— >oo N 

showing that the scattering length of N 2 V(Nx) is clq/N. Equivalently, this follows from 
/ dxN 2 V{Nx)f{Nx) = 8vra /iV. 

It follows immediately from Lemma 13.21 that < w(Nx) < c for some c < 1 and for all 
x£l 3 , and that there exists C with 

C N 

w{Nx) -n\xJTT and IVMNx)1 - c n^xYTi - (54) 

We will use the solution f(Nx) of the scaled zero-energy scattering equation to approx- 
imate the correlations among the particles, arising on the microscopic scale. It is however 
important to keep in mind that these correlations are also modulated on the macroscopic 
scale. The macroscopic variation is described, or at least, this is what we expect, by the 
solution of the modified Gross-Pitaevskii equation ()50p . We define therefore the kernel 

k t (x,y) = -Nw(N(x-y))<pl N \x)<pl N) (y) (55) 
and the corresponding unitary operator 

T(k t ) = exp Q J dxdy (k t (x,y)a* x a* -k t (x,y)a x a y ) 

In the next lemma, we collect several bounds for the kernel kt which will be useful in the 
following. 

Lemma 3.3. Let ip[ N ^ £ i^ 1 (1R 3 ) be the solution of pBp with initial data £ ii" 1 (M 3 ) . Let 
w(Nx) = 1 — f(Nx), where f solves the zero-energy scattering equation fl|). Let the kernel 
kt be defined as in A55\). 

(i) There exists a constant C, depending only on Wh 1 such that 

\\hh<c, 

||Vifct||2,||V 2 fct||2<CVW, 

\\V 1 (k t k t )\\ 2 ,\\V 2 (ktk t )h<C. 

Defining p(kt) and r(kt) as in (f^5[ ), so that ch{kt) = 1 + p(kt) and sh(kt) = kt +r(kt), 
it follows from Lemma \2.3l part (Hi) and (iv), that 

\\p(h)hAHk t )hA\sh(kt)h<c, 

\\Vip(kt)h,\\V2P(kt)h < C, 
\\Vir(kt)hA\V2r(h)h <C. 

(ii) For almost all x,y E M 3 , we have the pointwise bounds 

\h(x,y)\ <min (Ar|^)( x )||^)( y )|,^l^|^)( x )||^)( y) |^ , 

\p(h)(x,y)\ < Cl^^)!!^^)!. 
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(iii) Suppose further that <p £ H (R ). Then 

sup \\k t (.,x)\\ 2 , sup ||p(/c 4 )(.,x)|| 2 , sup \\r(k t )(.,x)\\ 2 , sup ||s/i(/c t )(., x)|| 2 < C\\ipf ] \\ H 2. 

xeK 3 xeM 3 xeM 3 



We will also need bounds on the time derivative of the kernels kt, p(kt), r(kt). These are 
collected in the following lemma. 

Lemma 3.4. Let ip G -fT 4 (lR 3 ), and ip[ € -ff 4 (R 3 ) be the solution of {5U\) . with initial data 
ip. Let w(Nx) = 1 — f(Nx), where f is the solution of the zero-energy scattering equation 
flj). Let the kernel kt be defined as in \55\). so that 

kt(x,y) = -Nw(N(x - y)) (4 N \x)<p { t N) (y) + <pf \x)$* \y)) . (56) 



Then there are constants C, K > 0, where C depends on the \\<p\\h a an d K only on ||</?||.ffi 
such that the following bounds hold: 

(i) 

Hfctlb, Hfeclb, ||p(fct)||2, 11^(^)112 < C7e« r [*l, 

(ii) 

l|Vip(A;)|| 2 , ||V 2 p(fe)|| 2 , ||Vir(A:)|| 2s ||V 2 f(A;)|| 2 < Ce*™, 

(iii) 

sup ||A; t (.,x)|| 2 ,sup ||p(fc t )(.,x)||2,sup||r(A;t)(.,x)||2,sup \\sh(k t )(., x)\\ 2 < Ce K ^. 

The proof of the last two lemmas can also be found in Appendix |Bj 

As explained in the introduction, we are going to approximate the many-body evolution 

e- mNt W(VN<p)T(k )i> 

of an initial state which is almost coherent, but with the correct short-scale structure, by 
the Fock state W{\fN(p\ N) )T(k t )^, which is again almost coherent and has again the correct 
microscopic correlations. This leads us to the fluctuation dynamics, defined as the two- 
parameter group of unitary transformations 

U(t;s) =T*{k t )W*{VNv[ N) y~ mN{t ~ s) W{VN<f { s N) )T{k s ) (57) 

where U(s; s) = 1 for all s £ R. 

The fluctuation dynamics satisfies the Schrodinger-type equation 

id t U(t;s) = C N {t)U{t;s) 

with the time-dependent generator 

C N (t) = T*(kt) [id t W*{VNip[ N) )\ W{jNip[ N) )T{k t ) 

+ T^k t )W%VN4 N) )n N W(VN<pl N) )T(k t ) + [id t T*{h)\ T(h). 

The next theorem, whose proof is deferred to Section [6j is the main technical ingredient 
of this paper. It contains important estimates for the generator C^i^t), which will be used in 
the next section to control the growth of the expectation of the number of particles operator 
with respect to the fluctuation dynamics U(t; s). 
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Theorem 3.5. Define the time- dependent constant (of order N) 
C N (t) = C { ^(t) + J dxdy\V x sh(k t )(y,x)\ 2 

+ J dxdy(N 3 V(N.)*\4 N) \ 2 ){x)\sh(k t )(y,x)\ 2 
+ J dxdydzN 3 V(N(x-y))i P ( t N \x)Tp ( t N \y)sh(k t )(z,x) sh(k t )(z,y) 
+ Re J dxdydz N 3 V(N(x - y))<p[ N) (x)<p { t N) (y)sh(k t )(z,x)ch(k t )(z,y) 
+ / dxdyN 2 V(N(x - y)) 



dz sh(kt)(z, x) ch(z, y) 



+ 



dz sh(k t )(z, x) sh(z, y) 



+ I dzidz 2 sh(k t )(z 1 ,x) sh(k t )(zi,y)sh(k t )(z2,x) sh(k t )(z 2 ,y) 



and let 

C N (t) = C N (t) - C N (t). 
Then we have, for some K > depending only on ||<£>||#i, 

,7V 2 



and 



Moreover, 



and 



c N (t) >\n N -c^- Ce K M (TV + 1) 



C N (t) <^<H N + C^- + Ce K ^ (Af + 1) 



± [Af, <U N + C^- + Ce K ^(N+l) 



±C N {t) <U N + Ce K M + 1 



4 Growth of fluctuations 



(58) 



(59) 

(60) 
(61) 
(62) 

(63) 



The goal of this section is to prove a bound, uniform in N, for the growth of the expectation 
of the number of particles operator with respect to the fluctuation dynamics. The properties 
of the generator £/v(£) of the fluctuation dynamics, as established in Theorem 13.51 pl a Y here 
a crucial role. 



Theorem 4.1. Suppose ip G T (possibly depending on N) with = 1 is such that 



(64) 



for a constant C > 0. Let ip € -£f 4 (R 3 ), and let ipi be the solution of the modified Gross- 
Pitaevskii equation |5Q|) with initial data if. LetU{t;s) be the fluctuation dynamics defined 



in (57). Then there exist constants C,c\,C2 > such that 

(ip,U*(t;0)MU(t;0)^} < Cexp(ci exp(c 2 |i|)). 
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The strategy to prove Theorem 14 . 1 1 consists in applying Gronwall's inequality. The deriva- 
tive of the expectation of M is given by the expectation of the commutator i[Af, where 
C]y(t) is the generator (fT4"j) of the fluctuation dynamics. By (|62|) . this commutator is bounded 
in terms of the energy, of (JV+ 1), and of M 2 /N (the difference between £jv(t) and the gener- 
ator >Cjy(i) appearing in (|62p is a constant and hence does not contribute to the commutator). 
The growth of the energy is controlled with the help of (|63p . What remains to be done in 
order to apply Gronwall's inequality is to bound the term M 2 /N. In the next proposition, 
we show that the expectation of M 2 /N at time t can be controlled by its expectation at 
time t = (a harmless constant, by the assumption (|64p ) and by the expectation of (A/" + 1) 
(which fits well in the scheme of Gronwall's inequality). 



Proposition 4.2. Let the fluctuation dynamics U(t;s) be defined as in (51). Then there 
exists a constant C > such that 

U*{t;0)Af 2 U{t;0) < C(NU*(t; 0)AfU(t; 0) + N(Af + 1) + (M + l) 2 ). 

The next lemma is useful in the proof of Proposition 14.21 

Lemma 4.3. Let k t £ L 2 (R 3 x R 3 ) be as defined in [55]) . Then there exists a constant C , 
depending only on ||^||2; such that 

T*(k t )MT(h)<C{M+l), (65) 

T*(k t )N 2 T{k t ) < C(N+1) 2 (66) 

for all t G R. 

Proof. We use the decomposition ch(fct) = 1 + p(kt) and the shorthand notation c x (z) = 
ch(fc()(z, x), p x (z) = p(k t )(z, x) and s x (z) = sh(k t )(z,x). We have 



(i>,T*(ktWT(h)*P)= I dx(ij,(a*(c x )+a(s x ))(a(c x ) + a*(s x ))ij) 

dx \\(a x + a(p x ) + a*(s x ))i)\\ 2 

< C j dx\\a x il>\\ 2 + j dx\\a{p x )il>\\ 2 + j dx\\a*(s 2 

< C(l + \\p(k t )\\ 2 + llsh^) 111)11^+ 1) 1/2 VH, 

and dSSD follows by LemmaES (since ||p(fct)|| 2j ||sh(Jfet)[| 2 < e} kt W' 2 ). To prove dnSJ) : we observe 
that 



{i>,T*(k t )N 2 T(k t )4>) = J dxdy (ij,T*(k t )a* x a x a* y a y T(k t )ij) 
dx{^,T*{k t )a* x Na x T{k t )i)) + (^,T* {k t )NT(k t )^) 

dx((a(c x ) + a*(s x ))^,T*(k)MT(k)(a(c x ) + a*(s x ))^) + ty, T*(k t )MT(k t )^). 

Then, applying ([65]) . we obtain 

{^,T*{k t )M 2 T{k t )^) 

<C j dx\\{M + l) 1/2 (a(c x ) + a*(s x ))4>\\ 2 + C(^, (TV + 

< C J dx {\\a x M l ' 2 n 2 + \\a(p x )N l ' 2 n 2 + h*(s x ){M + 2)V^f ) + C^, (Af + 

< C(l + HKMH 2 , + ||sh(A;i)|| 2 )^, (A/" + 1) 2 V). 

The bounds from Lemma 12.31 imply (|66p . □ 
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Proof of Proposition \4-S\ From Lemma 14.31 we nn d 

(ip,U*{t;0)N 2 U{t;0)^) < C(^,U*(t;0)T*(k t )M 2 T(k t )U{t-0)^). (67) 

We now show how to bound the r.h.s. of the last equation. Using the definition of the 
fluctuation dynamics W(t;0) = T*(kt)W*(VN(p^)e- m ^W*{y/Ntp)T(kQ), we find 

{iP,U*(t;0)T*(k t )N 2 T(k t )U(t;0)i{j) 

= {MT{k t )U(t;0)^,W*{^Nipi N) )(M - VNc/)^) + N)e- iHNt W(VNip)T(k )^) 

= {MT{k t )U{t^,W*{^ip ( ^)Me~ mNt W{VN^)T{k )^) (68) 

- y/N {NT{k t )U{t- 0)i/>, W*{y/N<pl N) ) 4){v[ N) )e- mNt W{VNv)T{kQ)i)) 
+ N{^M* (*; 0)T*(k t )MT(k t )U(t; 0)^) 

where we used the notation </>(/) = a(f) + a*(f), and the property ([36]) to show that 
W*(VNy\ N) )NW{VNy\ N) ) =M- jN<j){v\ N) ) + N. 

In the first term on the r.h.s. of (|68p . we use now the fact that N commutes with %n- in 
the second term, on the other hand, we move the factor cf)(ip[ N ^) back to the left of the Weyl 
operator W(^/N(p[ ), using that 

We conclude that 

(V#*(t; 0)T*(k t )M 2 T(k t )U(t; 0)^) 

= {MT(k t )U{t- 0)fj}, W*{Vn4 N) ) e~ mNt W(VNtp)(N + y/Nfttp) + iV) T(fc )^) 
- sfN {NT{k t )U{t-Q)i)j, (<K^ (A °) + 2^) W*(VN(p^)e- mKt W{VNip)T(ko)i/>) 
+ N(i;,U*(t;Q)T*(ktWT(k t )U(t;0)il>) 

= (NT(k t )U(t- 0)V>, W*(\ZiV>f°) e- iHNt W{VNip)NT(k Q )^) 

+ VN(NT(k t )U{t; 0)i/), W*{VN<p { t N) ) e- inNt W{^/N(p)^(<p) T{k )^) 
-y/N{M T{k t )U{t; 0)V>, <K¥?° )W*(VN<pi N) ) e~ W( v / iV^)T(A;o) V) . 

(69) 

By Cauchy-Schwarz, we obtain 

(ip,U*{t; 0)T*(k t )M 2 T(k t )U(t; 0)^) 

< \\MT{k t )U{t- 0)^|| 

x (||AAT(fco)V) || + v^V||<AM r(Ao)^|| + v^H^f ^TihMt- 0)V||) 

< |||ATT(fc t )ZY(i;0)^|| 2 

+ C (||AAT(A; )^|| 2 + iV||0M T(fc )^|| 2 + AT||<HrfV(^M*; 0)V>l| 2 ) • 

Subtracting the first term appearing on the r.h.s., and using the bound ||<£(/)?/>|| < ||/||2||(A/'+ 
1) 1/2 V>||, we find that 

(t; 0)T* (k t )Af 2 T(k t )U(t; 0)$) 

< c (pfTfaWf + N\\(Af + l) 1 /^)^ 2 + N\\(Af + l) 1 /^**)^ o)v>|| 2 ) . 
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By Lemma 14.31 we conclude that 

{i>M* 0)T* {k t )N 2 T{k t )U{t- 0)^> 

< CiV||AA 1/2 ^; 0)V|| 2 + C\\Af4>\\ 2 + CN\\(AT + l) 1/2 ^l| 2 - 

With (f67|) . this concludes the proof of the proposition. □ 

We are now ready to show Theorem 14. II 

Proof of Theorem \4-l\ Let C^it) be defined as in (|58[) and define 

U{t;s) = e^^^Wfts). 

Then U(t; s) satisfies the Schrodinger type equation 

id t U{t; s) = Z N (t)U{t; s), with U(s; s) = 1 

for all s E K, and with generator £jv(i) = £jv(i) — CV(t), as defined in ([59]) . On the other 
hand, since the two evolutions only differ by a phase, we have 

(ip,U*(t;0)NU(t;0)^} = {ip, U*(t; $)NU{t; 0)^). 

We now use the properties of £jv(i), as established in ([60]) . (|6T|) . ([62]) and ([63]) . Eq. (|60|) 
implies 

"Htv < 2£ 7V (t) + C^ + Ce^l'l(AA+l). (70) 

From Proposition l4.2l we conclude that there exist a constant C\ (depending on (ip, (N 2 /N + 
J\f+l)ip)), such that 

< (ip,U*{t-,0)H N U{t-,0)-ip\ < /ip,U*(t;0) (2£ N (t) + C ie m (N + l)^jU(t;0)^ . (71) 

From (|62p . combined with ()70p and Proposition 221 there exists moreover a constant C2 > 
(depending on (ip, (Af 2 /N +N + l)-0)) such that 

^ 0)/V"W(i; 0)^ < (ip,U*(t; 0) (aCtf(t) + C 2 e K ^(M + 1)) 0)V>) . 

We now estimate the growth of the expectation of the generator £jv(i). Using ([63]) together 
with (|7U|) and Proposition 14.21 we conclude that there exists a constant C3 > (again, 
depending on (ip, (Af 2 /N + M + 1)0}), with 

j t ^,W(t;0)C N (t)U(t;0)ip) = (iP,U*(t; 0)C N (t)U(t; 0)0 

< 0) (2£jv(t) + C 3 e K|4| (AA + 1)) 0)^) • 

We now fix D := max(Ci + 1, C 2 , C%,K). Then, we have 
j t ^,U* (t; 0) (£ N (t) + De*® (TV + 1)) 0)0 



< (iP,U*(t;0) [(2 + 2L>e^ |t| )£ JV (t) + (De x| * l + J D^e^ l ' l + J D 2 e 2 ^ l * l )(AA + l) W(t;O)0 
<(2 + J De^l'l) (iP,U*(t;0) (c N (t) + De K ^(AT + 1)) 0)0 

< 2Z)e x l*l <V, Z? (t; 0) (c N (t) + De*™ (/V + 1)) U(t; 0)V>) . 
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By Gronwall, we conclude that 

U,U*(t] 0) (c N (t) + De K ^(M + 1)) U(t;0)ip\ 
<Ce 2 # eK|t ' (V,(^v(0) + CAT+1))^) 

< Ce 2 ^ 141 ^, + + C{N + 1)) ^ 

where in the last inequality, we used the upper bound ([6T]) . From the assumption (}64"]) . we 
obtain 

U,U*(t;0) (c N (t) + DeFW(N + l))w(t;0)v) < Cexp(ci exp(c 2 |t|)). 
The claim now follows from (|70p . because D > C*x + 1. □ 

5 Proof of the main theorem 

Using the bounds established in Theorem 14. 1] we proceed now to prove our main result. 

Proof of TheoremEJl Let *jv,t = e- iHNt W{y/Nip)T(k )ip. The one-particle reduced den- 
sity of ^N,t has the integral kernel 

T N ] t( x ' y) = 7^ — \j^r-\ • ( 72 ) 

V*N,t,-NVN,t) 

We start by computing the denominator. Since commutes with the number of particles 
operator, we find 

(fi! N ,t,N*N,t) = (ip,T*(k )W*(VN<p)NW(VN<p)T(k )iP) 

= (V, T* (k ) (M - VN(f>(<p) + iv) r(fco)^) 

= JV + (V, r*(fc ) (N - VN<t>(<p)) T(k )^ . 

By Lemma 14.31 

(V, T* (k )AfT(k )ip) < Cty, (AT + 1) V) < C 

and 

|(V,T*(A;o)0( V )T(fc Q )V}| < C(V,r*(A:o)(AA+l) 1 / 2 T(fco)V) < C(V>, (AA + 1)V>> < C. 
Hence, there exists C > with 

|(^,AA^) - AT] ^CW 1 / 2 . (73) 

On the other hand, with ip[ denoting the solution of the modified Gross-Pitaevskii equation 
([50]) , the numerator of ([72]) can be written as 

= (if), T(k )W{VN^)e imN ala x e- mNt W{VN^)T(k )ifj) 
= N \p t {x)Tp t (y) 

+ VNip t (x) {TP,T*(k )W*(^Nip)e itHN (a* - Vn^\ N) (yj\ e- mNt W{^/N^)T(k Q )^) 
+ ^NTp t (y) {iP,T*(k Q )W*(VN(p)e itnN (a x - y/N^ (xj\ e- mNt W(VN<p)T(kv)ip) 
+ (y,T*(ko)W*(VN<p)e m » (a* y - VN^ t N) (y)) 

x (a x - VN<pl N) {x)) e~ inNt W{VNip)T(k )ijy 
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Recognising that 



(a* y - VNlpf\y)) = W(VN<pr'KW*(VN<pr>) 
(a x - ^N<pf\x)) = W(^/N<pl N) )a x W*(VN<pl N) ) 



we obtain 

{^N,t,a* y a x ^ Nt t) 



N^\x)lpf\y) 



+ y/N<p\ N \x) (i>,T*(k )W*{VN<py mN W(VN<p { t N >) 

x alW^VNip { t N) )e- mNt W(VNip)T(ko)^} 



+ VNTp ( t N) (y) (V, T* (k )W*{VNv)e mN W(VN^ t N) ) 

x a x W*(VNv { t N) )e- mNt W(VN V )T{k )ij) 



+ (^,T*(fc )W*(\/iVV)e iOTN W(^ 
Combining the last equation with (j73|) and inserting in (j72[) . we have that 

r%(x;y)-lpf\y)^\x) 
C 



< 



^ N) (x)U N \y)\ 



1 



+ Ho* W*(VNip { t N) )e- in ^W(VN^T(k 



+ 



4 N \x)\ \\ay W*(VN^>)e- l ' HNt W(VNcp)T(k 



'iV.Nt-1 



1 

iV 1 



+ — Way W* 



Nip { t N) )e- inNt W(VN<p)T(k 



x Wa x W*(VNp { t N) )e- mNt W(VNv)T(k 
Taking the square and integrating over x, y, we find 



Lemma 14.31 implies that 



< 



C 



W(M+1) 1/2 W*(VN(fl N) )e- iHNt W(VNtp)T(ko 



HS 



< 



c 



us 



N 
C 



W(N + 1) 1/2 r*(fc t )w*(\/F^ Ar) )e-^*w(\/iV^)r(M) 

wiM+i^u^mt 



with the fluctuation dynamics lA(t\ s) defined in (|57p . From Theorem 14. 1 1 we conclude that 



< 



Cexp(ci exp(c 2 |£|)) 



HS 



N 



Since l^^)^ 7 ^! is a rank-one projection, and T^L > 0, it follows that the difference 

— \(ft)((p t \ can have only one negative eigenvalue. Since the trace of — \(p t }{ip[ \ 
vanishes, it must have one negative eigenvalue, with absolute value equal to the sum of all 



27 



positive eigenvalues. As a consequence, the trace norm of the difference is controlled by the 
operator norm (given by the absolute value of the negative eigenvalue) and therefore also by 
the Hilbert-Schmidt norm. This shows that 



Tr 



< 



Cexp(ci exp(c 2 |i|)) 



Theorem 11.11 now follows because, if ipt denotes the solution of the Gross-Pitaevskii equation 
Proposition 13.11 implies that 



Tr 



i^i-i^rwf ! 



<2|| V 9 t -^ A ' ) || 2 <^^ 



□ 



6 Key bounds on the generator of the fluctuation dynamics 

In this section, we prove Theorem 13.51 concerning the generator £jv(i) of the fluctuation 
dynamics 

U(t;s) = T^k t )W*(VN^ N) )e- inN(t ~ s) W(VN^)T(k s ) 
as defined in (1571). We write 



(74) 



C N (t) = T*(k t ) id t W*(VN<p ( t N) ) W{y/N^ t N) )T{k t ) 
T*(k t )CP(t)T(k t ) + [id t T*(k t )]T(k 



+ T*(k t )W*(VN^ t N) )n N W(VNip { t N) )T(k t ) + [id t T*(k t )\ T(h) 



with 



c§\t) = \id t w*(VN<p { t N) )\ w(Vn<pI n) ) + w*(VN(fl N) )n N w(VN<p { t N) )- 

A simple computation shows that 

~id t W*(VN4 N) )] W(Vn4 N) ) = -a(VNid t 4 N) ) - a*(^Nid t ^ N) ) - N , id t ^ N) ) . 
On the other hand, using (|36p . we find 



w*(VNip { t N) )n N w(y/N(pl N) ) 



N 



\\Wt 



W||2 , 

2 ' 



a* ( (N*V(N.) * \4 N) \ 2 )^ N) )+a{ (N*V(N.) * 



J dxV x a* x V x a x + J dx(N 3 V(N.)*\^ t N) \ 2 )(x)a* x a x 

+ J dxdyN 3 V(N(x-y))^ t N) (x)Tp { t N \y)a* x a y 

+ \l dxdyN 3 V(N(x-y))(vl N \x)vl N \y)atal + lpl N) ^ 



+ 



+ 



dxdyN 3 V(N(x - y))a% ( rf%)a* + > (y)a y ) a a 



2N 



dxdyN 3 V(N(x — y))a x a*a y a x . 
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+ 



Combining the last two equations and using (|50p . we conclude that 
Z$(t) = N jdx (n*V(N.)(1 - f(N.)) * |^| 2 ) (x)U N \x)f 

+ VN [a* ((iV 3 «;(iV.)F(iV.) * \tp[ N) | 2 )^ (A °) + a ((N 3 w(N.)V(N.) * \tp[ N) l 2 )^ 
/" dxV x a* x V x a x + J dx (N 3 V(N.) * \(p[ N) \ 2 ){x)a* x a x 
+ J dxdyN3v{N{x-y))tp\ N \x)Tp[ N \y)a x a y 

+ \f d^iVV(iV( a; -y))(^ ) (^)^(y)a*a*+# ) (x)# ) (y)a :r a : 
+ J= y dxdy N 3 V(N(x - y))a* x ^ ] {y)a* y +Tp[ N) {y)a y ) a x 

where C^ N (t), for j = 0, ... ,4, is the part of (t) containing j creation and annihilation 
operators. Recall here that = 1 — /(x), as defined in ([52]) . 

From (|74p . we find that the generator of the fluctuation dynamics is given by 

4 

C N (t) = £® N (t) +Y J T*{h)C%{t)T(k t ) + [id t T*(k t )]T(k t ). (75) 
j=i 

In the next subsections, we study separately the different terms on the r.h.s. of (|75|) . The 
final goal of this analysis, a proof of Theorem 13.51 will be reached in Subsection 16.61 



Notation. In the rest of this section, we will use the shorthand notation 

c x (y) =ch(kt)(y,x), s x (y) = sh(k t )(y,x), p x (y) = p(k t )(y,x), r x (y) =r(k t )(y,x). (76) 

Moreover, ||p||2, IMb, [ | six ] 1 2 will denote the L 2 -norms of the kernels p(k t )(x,y), r(k t )(x,y), 
and sh(A^)(a;, y) over R 3 x R 3 (in other words, they denote the Hilbert-Schmidt norms of the 
corresponding operators). The norms \\px \ I 2 j 1 1 \ I 2 •> 1 1 sh^ 1 1 2 ; on the other hand, indicate norms 
over R 3 . Finally, the notation {.,.) will denote the I? inner product. We will abbreviate 
T{k t ) by T. 

6.1 Analysis of T*C [ ^ N {t)T 

Conjugating the linear term Cj^-At) produces again linear terms. From Lemma 12.31 we 
obtain 

T*cf ] N {t)T 



N J dxdyN 3 V(N(x-y))w(N(x-y))\^ N \y)\ 2 (^?\x)T* a* x T + Tp[ N \x)T* a x T 
N J dxdyN 3 V(N(x-y))w(N(x-y))\ip { t N \y)\ 2 ^ t N \x) (a*{c x ) + a(s x )) 
+ v / iV / dxdyN^ViNix-y^wiNix-y))^^ 2 ^) (a(c x ) + a* (s x )). 



(77) 
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These terms are potentially dangerous because they are large (of order v N) and do not com- 
mute with the number of particles. We will see however that they cancel with contributions 
arising from the cubic part C^ N (t). 

6.2 Analysis of T*cf^ N (t)T 
We write C { ^ N (t) =JC + C { ^ N (t), with 

K = j dx V x a* x V x a x 

being the kinetic energy operator, and we consider separately the effects of tC and of the 
other quadratic terms collected in C^ N {t). 

6.2.1 Properties of T*KT 
We have 

T*/CT = J dxV x (a*(c x ) + a(s x ))V x (a(c x ) + a*(s x )) 

= j dxV x a* (c x )V x a(c x ) + J dxV x a*(c x )V x a*(s x ) 

+ J 'dxVMs x )V x a(c x) + /W W V^ + J dx\\V x s x \\l 

Following dlSD , we decompose c x (y) = 5{x — y) + p x (y) and s x (y) = k(x,y) + r x (y). Hence 

T*KT = K, + f dxdy\V x sh kt (y,x)\ 2 

+ j dxV x ala(V xPx ) + j dxa*(V xPx )V x a x + j dxa* (V xPx )a(V xPx ) 

+ J dxV x a* x a* (V x k x ) + j dxa* (V xPx )a* (V x k x ) + J dxV x a* x a*(V x r x ) 

+ J dxa*(V xPx )a*(V x r x ) + j dx a(V x k x )V x a x + j dxa(V x r x )V x a x 

+ J dxa(V ' x k x )a(V ' xPx ) + j dxa(V x r x )a{V xPx ) + j dx a* (V x k x )a(V x k x ) 

+ j f dx a*(V x r x )a(V x k x ) + J dx a* (V x k x )a(V x r x ) + j dx a* (V x r x )a(V ' x r x ). 

(78) 

The properties of T*fCT are summarized in the next proposition. 
Proposition 6.1. We have 

T*KT = j dxdy\V x sh kt (y,x)\ 2 + fC 

+ W dxdy(Aw)(N(x-y)) U N) (x)^ t N) (y)a* x a* y + # ) (x)#° (y)a x a y ) ^ 



+ S K (t) 
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where the error Exii) is an operator such that for every 5 > there exists a constant Cs > 
with 



±£ K {t) < SIC + C 5 e m (TV + 1), 
± [TV, £ K (t)\ <5K, + C s e K ^ (TV + 1), 
±£ic(t) <5K + C h e m (TV + 1). 

To prove Proposition 16.11 we will use the next lemma. 

Lemma 6.2. Let j±,j 2 G L 2 (M 3 x R 3 ). Let ji x (z) '■= ji(z,x) for i = 1, 2. T/ien we /iaue 



(80) 



dx{'ip,a?(ji )X )a?(j2,x)'il>) 



(81) 



i/ere and m the following a" can 6e either the annihilation operator a or the creation operator 
a* . Moreover, for every 8 > 0, there exists C$ > such that 



dx{ip,V x ala i (ji tX )'ip) 
dx(ip,a i (j ltX )S/ x a x 'ip) 



< 5fC + C$11.71 1|| II (TV + 1) 1//2, 0|| 2 and, by conjugation 
<5/C + C 5 ||j 1 ||2||(TV+l) 1 /2 V; || 2 . 



(82) 



Terms where the argument of a creation and/ or annihilation operator is the kernel V x k x 
(whose L 2 -norm diverges as N — > oo) can be handled with the following bounds. For every 
5 > there exists Cs > s.t. 



dx(ip, a* (V x k x )J(jx tX )i(j) 
dx{ip,J(j ljX )a(V x k x )i(j) 



< 5K + C 5 {1 + |b'i|||)||(TV+ 1) 1/2 ^|| 2 and, by conjugation 
<<5AC + C (J (1 + |U 1 ||1)||(JV + 1) 1 /V|| 2 - 

(83) 



Moreover, we have 



dx{4>, a* (V ' x k x )a(V ' x k x )^) 



< C||TV 



(84) 



To control the time derivative of Exit), we will also use the following bounds. For every 
5 > there exists Cs > such that 



dx(ip,a*(V x h x )a i (ji, x )'4>) 



<5K + C 5 e K \\l + || J1 || 2 )||(/V+1) 1 



/ 2 ,/,ll 2 



Moreover, 



dx(i(},a*(V x k x )a(V x k x )ip) 
Proof. To prove (JHTJ) , we compute 



^Ce^HTV 1 / 2 ^!! 2 . 



(85) 



(86) 



dx(^a tt (ji iX )a tt (j2 i:c )V'> 



< / dx\\J(j liX )^\\ \\at(j 2 , x 

< (J dx\\j 1>x \\ 2 \\j 2 , x \\2 \ || (W + 1) 

^lljilbll^lblK/v + i) 1 / 2 ^!! 2 . 
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Eq. ()82p . on the other hand, follows by 



dx(4>,V x a* x J(j ltX )i>) < j dx\\V x a x iP\\\\a s (J 1}X )i>\\ 

<5(i>,ici>) + Cs J dxiu^iiiiKAA + i) 1 / 2 ^!! 2 
<^,x;v) + c l5 ||i 1 ||i||(Ar+i) 1 /2^||2. 

To show ([83]) . we need to integrate by parts. We write 

J dxa*(V x k x )J(j ltX ) = J dxdyV x k(y,x) a*J(ji }X ) 
and we observe that 

V x k(y,x) = -V y k(y,x) - Nw(N(y - x)) (Vip { t N) (x)^ (y) + ^ (x)W t N) (y) 
Hence 

J dxa*(V x k x )J(j ltX ) = J dxdyk(x,y)V y a* y J(j ltX ) 

dxdyNw(N(x - y))V<p { t N) {x)ip\ N) (y) a* y J(j 1>x ) 
dxdyNw(N(x - y))V<p { t N) (y)<p { t N) (x) a* y 
This implies, using Lemma 13.21 to bound Nw(N(x — y)), 
dx{ip,a*{V x k x )J{ji tX )%l}) 
< / dxdy\k(x,y)\\\V y a y ip\\ \\a^(j 1>x 



(87) 



+ C dxdy ■ 



< 



x-y\ 
5 J dxdy\4 N \x)\ 2 \\V 



\V<pW(x)\\4 N \y)\ + \v4 N \y)\U N) (x)\) \\a y ^\\ \\a*(h, x 



+ C J ^y-^I^^Plli^lllllKAT + ^^Vf + C J dxdy\V<p\ N \x)\ 2 \\a y n 2 
+ C j dxdyj^-^\4 N \x)\ 2 \\a y i>f + C j dxdy\V^ N \y)\ 2 \\j ltX \\ 2 \\(M + if'^f. 

Using Hardy's inequality, we conclude that for every 5 > there exists Cs > (depending 
on ||j'i||2, <5, llv 9 *^ Wh 1 ) such that 



(tp, / dxa*{V x k x )J(ii )X )i>) 



To show (EH), we write 



dx a* (V ' x k x )a(V ' x k x ) 



< 6(ip, W) + C 5 (l + || j! Hi) ty, (N + m. 



dxdyidy 2 V x k(yi,x)V x k(y 2 ,x) a* yi a y2 
dyidy 2 g(y 1 ,y 2 )a* yi a y2 . 



-I 
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We have 

dy\dy 2 g(yx , y 2 ) a* yi a y2 ip) 

< / dy 1 dy 2 \g(yi,y2)\\\ay 



< dy 1 dy 2 \g(y 1 ,y 2 )\ 2 J ij dy\dy 2 \\a yi i>\\ 2 \\a y2 i' M2 
= \\ g \\ 2 \\M l l 2 n 2 - 

Moreover, from the definition (|55p of the kernel kt and from the bounds of Lemma 13.21 we 
find 



1/2 / \ 1/2 



g\\ 2 < C / dy\dy 2 dx\dx 



2- 



ki - yi| 2 ki - 2/2| 2 k2 - yi| 2 k2 - 2/2 1 2 



Fi - yi\\xi - y 2 \\x 2 - yi\\x 2 - y 2 \ 
\^\ Xl )f\^\x 2 )\ 2 



< C / dy\dy 2 dx\dx 2 ' 



x\ - y\\ 2 \xi - y 2 \ 2 \x2 - yi\ 2 \x2 - 2/2I 2 



+ c 1 ^^^^^f^'^^ 



\xi - yir\x 2 - j/2 1 

<C f d Xl dx 2 - ^-—\^{x X )\ % \ipf>{X2)\ 2 

J \x\-x 2 \ z 
+ c( S n P ldy^^ N \y)\^\\^\\% 1 
< C 

for a constant C depending only on the // 1 -norm of <p[ . The last two bounds prove (|84j) . 

The inequalities ([85]) . (f86|) can be proven similarly to ([83]) and (|53]): this time, however, 
the bounds will contain the norm Wh 1 i which is bounded by Ce K ^, as proven in 

Proposition 13.11 □ 

Proof of Proposition \6.1l We prove the first bound in (|80p . To this end, we observe that 
Lemma f6.2l can be used to bound all factors on the r.h.s. of (|78p (using the uniform estimates 
for ||Vip(A^)||i2(R3 X K3-j, ||Virjfe t ||x,2(K3 X M.3) from Lemma [373]) . with two exceptions, given by 
the term 

dx V x a* a* ( k x ) 



and its hermitian conjugate. To control (|88p . we use that, from (|55p . 

V^ 4 (y,x) = -N 2 Vw(N(x _ ) (^f } (y) - Nw(N(x - y))Vpf° (x)<p[ N) (y). 

Hence 

dxV x a* x a*(V x k x ) = -N 2 I dxdyVw(N(x - y))^[ N) {x)^[ N) (y)V x a* x a 



y 

(89) 

N I dxdyw(N(x-y))Vtpf\x)4 N \y)V x ala y . 
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The last term can be written as 

f dxV x a* x a*(jx) (90) 



withj(y, x) = —Nw(N(x—y))'S/ip[ N \x)ip[ N \y). Combining Lemma l3.2l and Proposition 1331 
we find that j £ L 2 (M 3 x M 3 ), with uniformly bounded norm. Hence, Lemma 16.21 implies 
that, for every 5 > 0, there exists C§ > with 

dx (V, V x a* x a*(jM < m 1/2 n 2 + C 5 \\{M + if'Hf. (91) 

The first term on the r.h.s. of (|89|) . on the other hand, can be written as 

-N 2 J dxdyVw(N(x - y)) V ( f\x) i p ( f\y)V x ala* y 

= N 3 [ dxdy(Aw)(N(x-y))^ t N \x)ip { t N \y)a* x a* y 
2 / j„j..T7...(Krr„ ..\\r7,Ji*0f„\,J.*0, 



The first contribution on the r.h.s. of the last equation is large and appears explicitly on 
the r.h.s. of (|79p (it will cancel later, when combined with other terms arising from C^ N 

and C^ N ). The second term, on the other hand, is an error; integrating by parts, it can be 
expressed as 



N 2 I dxdyVw(N(x-y))Vip { t N) (x)^ t N) (y)alal 

= - N [ dxdyw(N{x - y))Vip { t N) (x)V^ { t N) (y)a* x a* 



N j dxdyw(N(x-y))v4 N) (x)^ t N) (y)a* x V y a; 

dxV(p ( t N) (x)a* x a*(Nw(N(x - .)V<p[ N) ) + f dyV y a*a*(j y ) 



with j(x,y) = —Nw(N(x — y))Vip[ N \x)ip[ N \y). The second term is bounded as in ([91 
The first term, on the other hand, is estimated by 

dxVtp { t N \x){a x ij,a*(Nw(N(x - .))V tp { t N) )ip) 

< Csup \\Nw(N(x - .))Vc^ A0 || 2 ||(A/' + 1) 1/2 ^|| 2 (92) 

X 

^ch^ii^ikaa + i)^ 



't 

|2 



and (JHH). 

Also the second bound in (|80p follows from Lemma 16.21 In fact, when one takes the 
commutator of M with the terms on the r.h.s. of (|78p one either finds zero (for all terms 
with one creation and one annihilation operators, which therefore preserve the number of 
particles), or one finds again the same terms (up to a possible change of sign). This follows 
because, by the canonical commutation relations 

[AT, o(/)] = -a(f) and [Af, a*(f)\ = a*(f) 
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for every / G L 2 (R 3 ). Finally, the third bound in (|80p is a consequence of Lemma 16.21 as 
well. In fact, the time derivative Sk^) is a sum of terms very similar to the terms appearing 
on the r.h.s. of (|78p . with the difference that one of the appearing kernels contains a time- 
derivative. Combining the estimates from Lemma [6.2l (including, in this case, also (|85p . (|86p ) 
with the bounds for ||ViPfc t ||2, ||Viffc t ||2 from Lemma 13.41 and with the bound for ||^ 
from Proposition 13.11 (needed to control terms similar to (|89p . (|92p . with a factor of iff 
replaced by <p\ ), we obtain the last inequality in (f80j) . □ 

6.2.2 Properties of T*Z^ N (t)T 

We consider now the other quadratic terms, collected in £2,aK*)> defined by C^ N (t) = 
JC + C^ N (t). We have 



T*C%(t)T= J dx{N^V{N.)^\^f) ( x )( a *(c x )+a(s x ))(a(c x ) + a*(s x )) 

+ f dxdyN 3 V(N(x-y))^ N \x)^ N \y)(a*(c x )+a(s x ))(a(c y ) + a*(s y )) 



+ i y dxdyiV 3 y(A^(x - y)) 

x ) (x)<pf^(y)(a*(c a .) + a(s a! ))(a*(c !/ ) + a(s y )) + h.c. 
Expanding the products, and bringing all terms to normal order, we find 
T-£%{t)T 

-/d,(TO).|^l')W 

x [a* (0^)0(00;) + a*(s x )a(s :c ) + a*(c x )a*(s :c ) + a(s x )a(c x ) + (s x , s x )] 
+ J dxdyN^iNix-y^i^ix^iy) 

x [a^c^a^) + a*(s y )a(s x ) + a*(c x )a*(s y ) + a(s x )a(c y ) + (s x , s^}] (93) 
+ dxdyN 3 V(N(x-y))^ t N \x)<p { t N \y) 

x [a*(ca.)a*(Cj,) + a*(c x )a(sy) + a*(c y )a(s x ) + a(s x )a(s y ) + (s x ,c^)] 
+ dxd y N*V{N{x-y))Tp\ N \x)Tp\ N \y) 

x [a(ca.)a(cj,) + a*(s y )a(c x ) + a*(s x )a(c y ) + a* (s x )a* (s y ) + (c y ,s x )] . 

The properties of T*cf ] N (t)T are summarized in the following proposition. 
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Proposition 6.3. We have 
T*Z^ ] N {t)T 



dx(N 3 V(N.) *\4 N) \ 2 )(x)(s x ,s x ) 
+ I dxdyN 3 V(N(x-y))^ t N \xp { t N \y)(s x ,s y ) 



(94) 

+ Re I dxdyN 3 V(N(x - y)) V [ N) (x) V [ N) (y)(s x , c y ) 



1 



+ ^ J dxdyN*V{N{x-y)) [^\x)^\y) a%a* y +# ) (x)# ) (y) a x a y 
+ S 2 {t) 

where the error 82(f) is such that, for appropriate constants C, K > 0, 

±S 2 (t) <Ce x l*l(AA+l), 
±[Af,£ 2 (t)] <Ce K ^(M+l), (95) 
±£ 2 (t)<Ce K ^ (Af+1). 

To show Proposition 16.31 we will make use of the next lemma. 

Lemma 6.4. Let ji,j 2 G L 2 (M. 3 x M 3 ). Let ji, x (z) := ji(z, x) for i = 1, 2. TTien i/iere exists 
a constant C such that 

^^^(^(x-y))!^^)!!^^!!!^^!^^!!^^^)^!! / n 

(96) 

^cil^ll^lbilbllislblKAA + i) 1 ^!! 2 . 

Moreover, 

dxo!y^V(Ar(x-y))|^ ) (x)||^ ) (y)|||a«(j 1 , :c )^||||a^|| 



(97) 
and 

dxdyN*V{N{x-y))\v?\x)\\v?\y)\\\a x m«vn < C\\vi N) \\ 2 H *\W +l) 1/2 i>f. (98) 



TTie bounds remain true if both creation and/ or annihilation operators act on the same 
variable, in the sense that 



dxdyN 3 V(N(x - y))\4 N) (x)\ 2 \\J(JiM\ hKh,y 

^cil^ll^lljilbll^lblKAA + i) 1 /^!! 2 , 

j ^rf2/JV 3 V-(7V(^ — 2/ ))|^ 7V) (^)| 2 ||a,ttO- 1 ,, / )^|| < C||^ 7V >||^ 2 ||j 1 || 2 ||(^-K l) 1 / 2 ^!! 2 , 

(99) 
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Proof. To prove (j96f) . we notice that for any a > 

cted^^CiNTCa: - 2/))!^ (a;)| |^(y)|[|a»C7i,x)^H l|a»(i 2 ,v) 
<a| ^ArM^-y))l^ ) (^)| 2 |U2, ? ;||ill(AA + l) 1 / 2 ^|| 2 

+ a' 1 1 dxdyN*V{N(x - y))\^(y)\ 2 \\jiM\^ + 1) V V|| 2 
< ||iVV(iV.) * irf^lU (a[|ji[|l + a-l^lll) ||(^+ lJ^Vf. 

Using ||iV 3 T/(iV.) * \tpt \ 2 \\oo < Cll^t ll^, and optimizing over a > 0, we obtain ([96]) . 
Analogously, (|97p follows from 

da:dyJVV(iV(x-y))|^(x)||^(y)|||a«OM)^IIK^II 

dsdyJVV^Ca-y^l^Ca;)! 2 !!^!! 2 

< ||ArM^-)*l^ JV) | 2 ||oo(a + a- 1 ||ji||l)||(JV'+l) 1 /Vf 

for any a > 0. Optimizing over a gives (|97p . Eq. (|98|) follows from Cauchy-Schwarz. The 
bounds in ()99p can be shown similarly. □ 

Proof of Proposition \6.3l To prove the first bound in (|95p , we notice that the quadratic terms 
on the r.h.s. of (j93|) can be controlled with Lemma 16.41 decomposing, if needed, a(c x ) = 
a x + a(p x ) and then applying ([96]) . ([97|) . or ([98]) . There are two exceptions, given by the 
terms proportional to a*(c x )a*(c y ) and its hermitian conjugate, proportional to a(c x )a(c y ). 
For these two terms the bounds from Lemma 16.41 do not apply. Instead, using a*(c x ) = 
a* + a*(p x ), we write 



/ 



dxdyN 3 V(N(x - y))^(z)^(vKfeK(cy) 

dxdyN 3 V(N(x - y))<p[ N \ X )<p[ N \y)a* x a* y 
+ j dxdyN 3 V(N(x-y))^ t N \x)^ t N) (y)a*(p x )a* y 
+ / dxdyN 3 V(N(x-y))^ t N \x)<p { t N) (y)a*(c x )a*( Py 



(100) 



The contribution of the last two terms can be bounded by Lemma 16. 4[ because one of the 
arguments of the creation operators is square integrable. In fact 

dxdyN 3 V(N(x - y))v[ N) (x)^ (y)(V, a* {pM) 



< J dxdyN 3 V(N(x-y))\^ N \x)\\^ N) (y)\\\a y ^\\\\a*(p x 



/2„/.i|2 



by ()97p and ()5ip . and similarly for the last term on the r.h.s. of (jlOOp . The hermitian 
conjugate of (|100|) . proportional to a(c x )a(c y ), can be handled identically. 
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The second bound in (|95p follows similarly, using the fact that the commutator of N 
with the terms on the r.h.s. of (|93p leaves their form unchanged (apart from the constant 
terms and the quadratic terms with one creation and one annihilation operators, whose 
contribution to the commutator [M, £2^)] vanishes). 

Also the third bound in (|95p can be proven analogously, using the bounds for ||shfc t ||2 
and llpfcjb) a s proven in Lemma 13.41 , When the time derivative hits the factor (p\ N \x) or 
<Pt (y), h generates a contribution which is bounded by \\tpt \\h 2 ^ Cllv 3 ^ ||1L < Ce K ^ 
(for some K depending only on here we used Proposition 13. ip . □ 

6.3 Analysis of T*C { ^ N {t)T 

We consider now the contributions arising from the cubic terms in Cx N (t). We have 

T*cf ] N {t)T 
1 



dxdyN 3 V(N(x - y)) 

tp[ N \y) (a*(c x ) + a(s x ))(a*(c y ) + a(s y ))(a(c x ) + a*(s x )) + h.c. 
)= j dxdyN 3 V(N(x-y))^ t N \y) 



x [a*(c x )a*(c y )a*(s x ) + a* (c x )a* (c y )a(c x ) + a* (c x )a(s y )a* (s x ) + a* (c x )a(s y )a(c x 
+a(s x )a*(c y )a*(s x ) + a(s x )a* (c y )a(c x ) + a(s x )a(s y )a* (s x ) + a(s x )a(s y )a(c x )] 

+ h.c. 

Writing the terms in normal order, we find 
T*A{t)T 



-j= J dxdyN 3 V{N{x-y)) V f\y) 

x [a* (c x )a* (c y )a* (s x ) + a* (c x )a* (c y )a(c x ) + a* (c x )a* (s x )a(s y ) + a* (c x )a(s y )a(c x ) 
+a*(c y )a*(s x )a(s x ) + a* (c y )a(s x )a(c x ) + a*(s x )a(s x )a(s y ) + a(s x )a(s y )a(c x )] 

+ 4= J dxdyN 3 V(N(x - y))^ N \y) 



x [{s y , s x )(a*(c x ) + a(s x )) + (s x , c y ) (a(c x ) + a*(s x )) + (s x , s x ) (a*(c y ) + a(s y ))] 
+ h.c. 

(101) 

The properties of T*C^ N (t)T are summarized in the following proposition. 
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Proposition 6.5. We have 



ryi* /* (0) rji 

1 's^ 1 



1 



,^ j dxdyN 3 V{N(x - y)) 

x <Pt(y)kt{x,y) {a(c x ) + a*(s x )) +Tpf\y)k t {x,y) (a*(c x ) + a(s x )) 
+ S 3 (t) 

-VN f dxdyN 3 V(N(x - y))w(N(x - y))\^ N) {y)\ 2 ^\x){a{c x ) + a*{s x )) + h.c. 



+ £ s (t) 

(102) 

where we used the definition $55)) of the kernel kt and where the error term £%(£) is such that 
for every 5 > there exists a constant Cs > with 



±£ 3 (i) <5 J dxdyN 2 V(N(x - y))a* x a* y a x a y + $ + C S e m (M + 1) , 

± W, £s(t)} < 5 J dxdyN 2 V(N(x - y))a* x a* y a x a y + + C s e K ^ (Af + 1) , (103) 

f N 2 
±4(0 < $ / dxdyN 2 V{N{x - y))a* x a* y a x a y + $ + C 5 e m (M + 1) . 



Notice here that the first term on the r.h.s. of (|102j) cancels exactly with the contribution 
(|77|) ; we will make use of this crucial observation in the proof of Theorem 13.51 below. 

Proof. To bound the cubic terms on the r.h.s. of (|10ip . we systematically apply Cauchy- 
Schwarz. This way, we control cubic terms by quartic and quadratic contributions, which 
are then estimated making use of Lemma 16.41 fthe quadratic part) and Lemmas 16.71 and 16.81 
(the quartic part). For example, 
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dxdyN 3 V(N(x - y))^ t N) (y)(ij,a*(c x )a*(c y )a*(s x )4>} 

<-Ly dxdyN 3 V(N(x - y))\ipi N \y)\\\a(c x )a(c y )i;\\\\a*(s a 



<^J dxdy N 3 V(N(x - y))\\a(c x )a(c y 

+ C S J dxdyN 3 V(N(x-y))\^ t N \y)\ 2 \\a*(s x 

<™J dxdyN 3 V{N(x-y))\\a y aM 2 + C 5 W + l) l ^\\ 2 

where, in the last line, we used (|111|) (from Lemma |6.8|) and (|99|) (from Lemma |6.4|) . All 
other cubic terms can be bounded similarly. We always separate the three creation and/or 
annihilation operators putting a small weight 5 in front of the quartic term and in such a 
way that, in the resulting quartic contribution, two operators depend on the x and two on 
the y variable. The corresponding quadratic term depends on x and can always be bounded 
by (f99|) . It should be noted that the quartic contribution has either the form ||a(c a; )a(c 2y )V , || 2 
or ||a(c x .)a"(j ?/ )V'|| 2 or \\a^(ji :X )a^(j2, y )ip\\ 2 , with square-integrable arguments ji,j2 (here a" 
is either a or a* ) . These terms can always be controlled using Lemma 16.71 or Lemma 16.81 As 
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for the linear contributions on the r.h.s. of (jlOip . the first and third can simply be bounded 
by N~ 1 ' 2 {M + l) 1 / 2 , since 

\(s y ,s x )\<C\^ N \x)\\^ t N \y)\. 
To bound the second linear term, we write 



1 



dxdyN 3 V(N(x - y))^ t N) (y)(s x ,c y )(a(c x ) + a*(s x )) 

j= J dxdyN 3 V(N(x-y))4 N) (y)kt(x,y) (a(c x ) + a* (s x )) + £(t) 



where ±£(t) < N-^ 2 (Af + 1) 1/2 because, using Lemma 13.31 

\(s x ,cy) -k t (x,y)\<C\vl N \x)\\<pl N \y)\. 

From (|5ip , this concludes the proof of the first estimate in (|103j) . The other two estimates 
are proven analogously, using the fact that the commutators of JV with the terms on the 
r.h.s. of (jlOip have the same form as the terms on the r.h.s. of (|10ip (with possibly just a 
different sign), and using the bounds for ||shfc t ||2 and HpftJIa from Lemma [374~|. and the bound 
for \\ip[ N ^ \\jj2 from Proposition 13.11 □ 

6.4 Analysis of T*C^ N T 

We consider next the contributions arising from the quartic part C^ N of C^\t). We have 

T*C%{t)T~ j dxdyN 3 V(N(x-y)) 

x (a*(c x ) + a(s x ))(a*(cy) + a(s y ))(a(c y ) + a* (s y ))(a(c x ) + a*(s x )). 
Expanding the products, we find 
2T*£,fh(t)T 

dxdy N 2 V(N(x - y)) a* (c x )a* (c y )a* (s y )a* (s x ) + a* (c x )a* (c y )a* (s y )a(c x ) 

+ a* (c x )a* (c y )a(c y )a* (s x ) + a* (c x )a* (c y )a(c y )a(c x ) + a* (c x )a(s y )a* (s y )a*(s x ) 
+ a* (c x )a(s y )a* (s y )a(c x ) + a* {c x )a(s y )a(c y )a* (s x ) + a* (c x )a(s y )a(c y )a(c x ) 
+ a(s x )a* (c y )a*(s y )a*(s x ) + a(s x )a* (c y )a* (s y )a(c x ) + a(s x )a* (c y )a(c y )a* (s x ) 
+ a(s x )a* (c y )a(c y )a(c x ) + a(s x )a(s y )a* (s y )a* (s x ) + a(s x )a(s y )a* (s y )a(c x ) 

+ a(s x )a(s y )a(cy)a* (s x ) + a(s x )a(s y )a(c y )a(c x ) . 
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Writing all terms in normal order, we obtain 
2T*C%(t)T 



(104) 



dxdy N 2 V(N(x - y)) a* (c x )a* (c y )a* (s y )a* (s x ) + a* (c x )a* (c y )a* (s y )a(c x ) 

+ a* (c x )a* (c y )a* (s x )a(c y ) + a* (c x )a* (c y )a(c y )a(c x ) + a* (c x )a* (s y )a* (s x )a(s y ) 
+ a* (c x )a* (s y )a(s y )a(c x ) + a* (c x )a* (s x )a(s y )a(c y ) + a* (c x )a{s y )a(c y )a(c x ) 
+ a* (c y )a* (s y )a* {s x )a(s x ) + a* (c y )a* (s y )a(s x )a(c x ) + a* (c y )a* (s x )a(s x )a(c y ) 
+ a* (c y )a(s x )a(c y )a(c x ) + a* (s y )a* (s x )a(s x )a(s y ) + a* (s y )a(s x )a(s y )a(c x ) 

+ a* (s x )a(s x )a(s y )a(c y ) + a(s x )a(s y )a(c y )a(c x ) 

+ J dxdy N 2 V (N (x - y)) (c y , s x )a* (c x )a* (c y ) + (s x ,c y )a(c y )a(c x ) 

+ 2( s y , s y )a* (c x )a* (s x ) + 2(s y , s y )a(s x )a(c x ) + 2(s y , s x )a* ic x )a* (s y ) 
+ 2(s x ,s y )a(s y )a(c x ) + 2{s y , s y )a* (c x )a(c x ) + 2(s y , s x )a* (c x )a(c y ) 
+ 2( Sxj s y )a* (s x )a(s y ) + 2(5^, s y )a* (s^a^) + (cy, s x )a* (c x )a(s y ) 

)a*(s y )a*(s x ) + ( Cj/j s x }a(s x )a(s y ) 

+ (s x ,c y )a*(s x )a(c y ) + (c y , s x )a*(c y )a(s x ) 
dxdyN 2 V(N(x - y)) [\ (s x , c y )\ 2 + \{s x ,s y )\ 2 + (s y , s y )(s x , s x )] . 

The properties of T*C^ N T are summarized in the next proposition. 
Proposition 6.6. We have 

2T*C^ N (t)T = J dxdyN 2 V{N{x - y)) [\(s x ,c y )\ 2 + \(s x , s y )\ 2 + (s y , s y )(s x , s x )] 

+ J dxdyN 2 V(N(x - y))a* x a*a y a x (105) 

+ / dxdy N 2 V(N(x - y)) (k(x, y)a* x a* + k(x, y)a x a y ) + £±(t) 



where the error £/±{t) is such that, for every 5 > 0, there exists a constant Cs > with 

M 2 



±£i(t) < 5 J dxdyN 2 V(N(x - y))a* x a* y a y a x + C s — + C s e*W (JV + 1) , 

f N 2 
± [TV, £ A (t)\ <5 dxdyN 2 V{N(x - y))a* x a* y a y a x + C 5 — + C s e m (M + 1) , (106) 

±£i(t) <S J dxdyN 2 V{N(x - y))a* x a* y a y a x + C s e m +M+ l^j . 

To prove Proposition 16. 6( we will make use of the following two lemmas. 
Lemma 6.7. Suppose j\,j2 £ L 2 (K 3 x M 3 ) are kernels with the property that 



M; := max sup / dy \ ji(x, y)\ 2 , sup / dx \ ji(x, y)\ 2 ) < oo 
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for i = 1,2. Let ji, x (z) '■= ji(z,x) and recall the definition c x (z) = chk t (z,x) from f"75| ). 
Then there exists a constant C depending only on Mi,M 2 and on the L 2 -norms \\j1W2, HJ2II2 
such that 

J dxdyN 3 V{N{x - y))\\J(j ljX )a y ipf < CMi\\(M + l)tp\\ 2 (107) 

and 



dxdyN 3 V(N(x - y))\\J(j 1>x )a}(j 2 , y W < C (Mx\\j 2 \\l M 2 \\j x \\l) \\ (M + 1)^|| 2 . 

(108) 

As a consequence 

dxdyN 3 V(N(x - y))\\a}(j^ x )a{c y )il,\\ 2 < CM 1 \\(M + 1)V>|| 2 . (109) 

The inequalities remain true (and are easier to prove) if both operators act on the same 
variable. In other words 



dxdyN 3 V(N(x - y))\\a}(j ltX )a x ^\\ 2 < CM 1 \\(M + 1)V>|| 2 , 
J dxdyN 3 V(N(x - y))\\J(j hx )a\j 2 , x W < C ^\n{M l \\j 2 \\ 2 ,M 2 \\j 1 \\ 2 )\\{M + I 
dxdyN 3 V(N{x - y))\\aHji, x )a(c y )^\\ 2 < CM^M + 1)^|| 2 . 



(110) 



Here 0} is either the annihilation operator a or the creation operator a* . 
Proof. To prove (|107p . we observe that 

dxdyN 3 V(N(x - y))||a s (i liX )a^|| 2 

< j dxdyN 3 V(N{x-y))\\h, x \\l\\{N+l) l ' 2 a y 

<Mi i dxdyN 3 V(N(x - y))\\a y M 1/2 ij\\ 2 



= CM 1 \\NiP\\ 2 . 

As for ()108p . we notice that (considering for example the case a^(j 2ty ) = a*(j 2jy )) 
dxdyN 3 V(N{x - yMJij^yih, 



< / dxdyN 3 V{N{x-y))\\hM\{N +l) 1/2 a*(j 2 



2 



< J dxdyN 3 V(N{x-y))\\j l X\\°*{h,yW + 2) 1 l 2 i,\\ 2 

< J dx^M^-y))IIJMllilU2,,ill(AA+l) 1 / 2 (AA + 2) 1 /V|| 2 

^cM.whwiuM+mi 2 . 

Eq. (|109p follows from the first two, by writing a(c y ) = a y + o,{p y ) (recall here that we 
are using the notation p y {z) = p(kt)(z,y) with the kernel p(kt) £ L 2 (M 3 x M 3 ) defined in 
Lemma f3. 3p , Eq. (jllOp follows similarly; in this case, however, one can immediately integrate 
over the variable y, simplifying the proof. □ 
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Terms of the form (|109p . but with ji x replaced by c x (which is not in L 2 ) are treated 
differently. 



Lemma 6.8. Recall the definition c x (z) = chk t (z,x) from ( 76). Then there exists a constant 
C > with 



dxdyN 3 V(N(x - y))||a( C;c )a( C2 ,)^|| 2 < C j dxdyN 3 V(N(x - y))\\a x a y ^\\ 2 

+ c\\(N + i)n 2 - 

More precisely, we have 

J dxdyN 3 V(N(x - y))\\a(c x )a(c y )ip\\ 2 = j dxdyN 3 V(N(x - y))\\a x a y ^\\ 2 +S(t) (112) 

where the error £{t) is such that, for every 5 > 0, there exists a constant Cs with 
±£{t) < 5 I dxdyN 3 V(N(x - y))\\a x a y ^\\ 2 + C 5 \\{M + l)-'-"' 2 



Proof. We write a(c x ) = a x + a(p x ), using the notation p x {z) = p(kt)(z,x) introduced in 
761). We have 



\\a(c x )a(cy)2p\\ < \\a x aytl)\\ + ||a a; o(p 2/ )V?|| + \\a(p x )a(cy 

Therefore, using (|107[) and ()109p . we immediately find (using Lemma 1331 to bound ||p||2 and 
sup x ||p*|| 2 ) 

J dxdyN 3 V{N{x-y))\\a(c x )a{c y )^\\ 2 <C J dxdyN 3 V{N{x-y))\\a x a y ^\\ 2 + C\\{N + 1 
To prove ()112p . we notice that 

dxdyN 3 V(N(x - y))\\a(c x )a(c y ^ 2 
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dxdyN 3 V(N{x - y))(ip, a* (c x )a* (c y )a(c y )a(c x )^} 
dxdyN 3 V(N(x — y))(t(i, a*a*a y a x ifj) 
+ J dxdyN 3 V(N(x - y))(ip, [a*(p x )a*a y a x + a* (c x )a* (p y )a y a x 

+ a* (c x )a* {c y )a(p y )a x + a* (c x )a* : (cy)a(cy)a(p x )]ij)^ 
=: j dxdyN 3 V{N{x - y))\\a x a y ^\\ 2 +£(*) 

where 

\£{t)\ < J dxdyN 3 V(N(x-y)) \\a(p x )a y ip\\\\a y a x 'ip\\ + ||a(c a; )a(pj / )'^'||||aj,a a .^[| 

+ \\a(c x )a(c y )ip\\\\a(p y )a x ip\\ + \\a{c x )a{c y )tp\\ \\a(c y )a(p XJ 
< 5 J dxdyN 3 V(N(x - y)) [\\a x a^\\ 2 + ||a( C;c )a( C?; )V|| 2 ] 

+ C 5 J dxdyN 3 V(N(x - y)) [\\a( Px )a y ^\\ 2 + \\a(c x )a(p y )^\\ 2 ) 

<5 j dxdyN 3 V{N{x - y))\\a x a y il>\\ 2 + C S \\{N + 1)V|| 2 - 
Here, in the last inequality, we used (|lU7p . (|1U9|) from Lemma I5T71 and (jllip . □ 
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Proof of Proposition \6.6\ To prove the first bound in f)106|) we observe that all quartic terms 
on the r.h.s. of ()104|) can be bounded using Lemmas 16.71 and 16.81 For example, the contri- 
bution arising from the first term on the r.h.s. of (|104p is bounded by 

dxdy N 2 V(N{x -y))(ij,a*(c x )a*(c y )a*{sy)a*{s x )ip) 

< J dxdyN 2 V(N(x - y))\\a(c x )a(c y )^\\\\a*(s y )a*(s, 

<5 J dxdyN 2 V{N{x-y))\\a{c x )a{c y )^\\ 2 

+ C 5 [ dxdyN 2 V(N(x - y))\\a*(s y )a*(s x ) 



<C5j dxdyN 2 V(N(x - y))\\a x a y i>f + ^-\\(Af + P-" 2 

where, in the last inequality, we used (jllip and (|108j) . All the other quartic terms on the 
r.h.s. of (|104p . with the exception of the fourth term (the one containing only c x or c y as 
arguments of the creation and annihilation operators), can be bounded similarly; the key 
observation here is that all these terms have at least one creation or annihilation operator 
with square integrable argument (this allow us to apply Lemma l6.7p . Moreover, in all these 
terms, the quartic expression does not contain the annihilation operators a(c x ) and a(c y ) in 
the two factors on the left, nor the creation operators a*(c x ) and a*(c y ) in the two factors 
on the right (in Lemma 16.71 in particular in (|1U9|) it is of course important that the factor 
a(cy) in the norm appears as an annihilation and not as a creation operator). To bound the 
fourth term on the r.h.s. of (|104p . where all the arguments of the creation and annihilation 
operators are not integrable, we cannot apply Lemma [6. 71 Instead, we use (|112p from Lemma 
6.81 We obtain 

dxdyN 2 V(N(x - y))(ip,a*(c x )a*(cy)a(cy)a(c x )tp) 

dxdyN 2 V(N(x - y))(if>, a* x a* y a y a x ^) + £(t) 

where the error £(t) is such that, for every 5 > 0, there exists Cg > with 

\£(t)\<6 [ dxdyN 2 V(N(x-y))\\a x a y n 2 + %\{N + l)n 2 - 



The quadratic terms on the r.h.s. of (|104j) can be bounded using Lemma 16.41 To this 
end, we observe that 

\(s x> s y )l\(s x ,c y )\<CN\d N \x)\\vl N \y)\. 

It is therefore easy to check that all the quadratic terms, with the exception of the first two 
(the quadratic terms appearing on the eighth line of (|1U4|) ). have a form suitable to apply 
one of the bounds in Lemma 16.41 More precisely, we apply (|96p . if the arguments of the two 
creation and/or annihilation operators are either s x or s y . If, on the other hand, one of the 
two arguments is c x or c y and the other one is s x or s y , we write a^(c x ) = a x + a*(jp x ) and 
then we apply ([96]) (to bound the contribution proportional to ar{p x )) and ([97]) (to bound 
the contribution proportional to a x ). Finally, if both arguments are either c x or c y (and we 
have exactly one creation and one annihilation operators), we write c$(c x ) = a x + cfi(p x ) 
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and we apply (|96|) . (|97|) and (f98|) . To control the two remaining quadratic contributions, we 
observe that, writing a*(c x ) = a* + a*(Pa;)) 



dxdyN 2 V(N{x - y)){c y , s x ){if>, a* {c x )a* (cy)rp) 

dxdyN 2 V(N(x-y))(c y ,s x )(i> 
+ [ dxdyN 2 V(N(x - y)){c y ,s x )(^,a*(p x )a* y ^) 

dxdyN 2 V(N(x - y))(c y ,s x )(ip,a*(c x )a*(p y )ij). 



(113) 



Since \(c y ,s x )\ < C N\<Pt (x)\\tpt (y)\, the last two terms can be bounded (in absolute 
value) using (f97|) and (f98|h respectively. We find 

dxdyN 2 V(N(x - y))(c y ,s x )(iP,a*(p x )a;iP) < C||^ (A ° ||^ 2 ||(Af + lf^f, 
dxdyN 2 V(N(x-y)){c y ,s x ){^,a*(c x )a*{p y )^) < C\\<p { t N) f H , \\ (M + 1) 1/2 VI| 2 - 

As for the first term on the r.h.s. of (|113p . we notice that 

(c y ,s x ) = k(x,y) +g{x,y) 
where g(x, y) = r(x, y) + {p y , s x ) is such that 

\g(x,y)\<C\^ t N \x)\\4 N) (y)\. 

Therefore, 

dxdyN 2 V(N(x-y))(cy,8 x )(1>,a*a*il>) = J dxdyN 2 V(N(x - y))k(x,y)(^,a* x a* y ^) 

+ J dxdyN 2 V(N(x - y))g(x,y)(^,a* x a* y ij). 

The second term can be bounded by 

dxdyN 2 V{N(x - y))g(x,y)(^,a x a*>ip) 



<C j dxdyN 2 V(N(x-y))\^ N \x)\\v[ N \y)\\\a x a y ^ 
<6 [ dxdyN 2 V{N{x-y))\\a x a y n 2 + C & [ dxdyN 2 V{N{x - y))\\^\x)\ 2 \^\y)f 



<5 j dxdyN 2 V(N{x - y))\\a x a y ?P\\ 2 + C 5 \\<p { t N) \\ 2 H i- 

Proceeding analogously to control the second quadratic term on the r.h.s. of (|104p . we con- 
clude that 

dxdy N 2 V(N(x - y)) [(c y , s x ) a* (c x )a* (c y ) + (s x , c y ) a{c y )a(c x )] 

dxdy N 2 V(N(x — y)) [k(x, y)a x a* + k(x, y)a x a y ] 
+ £(t) 
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where the error £ (t) is such that, for every 5 > there exists Cs with 

±S(t) <S J dxdyN 2 V(N(x - y))a* x a* y a y a x + C s \\y\ N) f H2 {M+l). 

We then use (|5"Tj) to conclude the proof of the first bound in (|106|) , 

The proof of the second inequality in (|106p is analogous, because commuting the terms 
contributing to £&(t) with the number of particles operator J\f either gives zero or leaves 
the terms essentially invariant (up to a constant and a possible sign change). Finally, also 
the third estimate in (|1Q6|) can be proven similarly, because the time derivative of the terms 
contributing to £^{t) can be expressed as linear combination of terms having the same form, 
just with one argument c x , Cy, s x or Sy replaced by its time-derivative. These terms can then 
be handled as above, using however the bounds for ||shfc t ||2 and ||pfc t ||2 from Lemma and 
the bound for H^^Hh 2 from Proposition 13.11 □ 

6.5 Analysis of [id t T*]T 
We set 

B = i j dxdy [kt(x,y)a%a* y - k t {x,y)a x a y ^ 



and 



B = \ dxdy (k t (x,y)a* x a y - k t (x,y)a x a y ^ 



with 

h(x,y) = -Nw(N(x-y))<pf\ X )vl N) (y) 

and 

k t (x,y) = -Nw(N(x - y)) (^ (JV) (x)^ ) (y) + ^{ N) {x)^{ N) {y) 
Then T = exp(B). Using (@7D, we find 

(d t T*)T= f d\e- XB VB(t)e XB ® = V--1t:4(B) (114) 
Jo ^o (n + 1)! 

up to an error term, which however converges to zero in expectation on the domain D(N) of 
the number of particles operator (this can be shown as in Lemma |2.3|) . Notice that, by the 
estimates in Prop. 16.101 the series is absolutely convergent in expectations. Since D(N) is 
invariant w. r. t. the fluctuation dynamics U(t,s) (this is proven similarly to Prop. W7Z\i , we 
can use (|114p to compute the expectation of (dfT*)T in the state U(t; 0)tp for any tp £ D(M). 
Next, we compute the terms on the r.h.s. of (|114p . 

2/"HI>3 v T®3\ 



(115) 



Lemma 6.9. For each n G N there exist f n ,i, fn,2 £ L (Ir x Mr) such that 

ad%(B) = i j dxdy (f njl (x,y)a*a* + f n ^(x,y)a x a y ) for all even n and 
adg(I?) = i J dxdy (f n ^{x,y)a* x a y + f n)2 (x,y)a x a*) for all odd n 

where 

\\fnA\2<2 n \\k t \\^\\k t \\ 2 , (116) 

for all n > and i = 1, 2, 

\ \\kt\\2 if n = 

\\fn,ih<{ An n, nn-1 f«t II lit L\\2\ ( 117 ) 



± n \\h\\T IIM2IIM2 + INI2 *f n > 1 
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and 

J dx\f nii (x,x)\ < 2lfc t ||2||fc t [| 2j J dx\f n!i (x,x)\ ^^Hfctlir 1 (\\h\\i + \\kh\\hh) 

(118) 

/or all n > 1. 

Proof. The proof is by induction in n. For n = 0, 

Hence fo,i(x,y) = k t (x,y) and fo,2(x,y) = k t (x,y), and the estimates (|1 16j) and (|117p are 
clearly satisfied. Suppose now the statement holds for some n S N. We prove them for 
(n + 1). We assume first that n is even. Then, using the canonical commutation relations 
we find 

ad« +1 (£) 

= [S,ad£(B)] 

^ y dxdy (k t (x,y)a* x al - k t (x,y)a x ay) ,^ f dxdy (f nil (x,y)a* x a* y + f n ,2(x,y)a x a y ) 

__ 1 
~ 2 

where 

(a:, 2) = ~\ j d V ( k t{x,y) (fn,2(z,y) + f n ,2(y,z)) + h(y,z) (/„, 2 (x,y) + /„, 2 (y, x))) > 

f n+li2 (x,z) = J dy (h(y,z) ( fn,i(x,y) + f n ,i(y,x)) + k t (x,y) (f n ,l(z,y) + /n,l(y> z ))) • 

(119) 



2 . 

dxdz (f n +i,i(x, z)a* x a z + / n+ i i2 (x, z)a x a* z ) , 



By Cauchy-Schwarz (similarly to (|49p ). we have 



<2||fc t || 2 ||/ n , 2 || 2 <2" +1 ^-» n+1 



||;n+l,l||2 ^ ^||«f||2||yn,2||2 ^ * ||«t|| 2 \\ K t\\2, 
||/n+l,2||2< 2||fc t || 2 ||/ n , 1 || 2 <2" +1 ||^||™ +1 ||^|| 2 , 

where we used the induction assumption. Moreover, again by Cauchy-Schwarz, 



(120) 



dx\f n+1 ,i(x,x)\ < 2||fc t || 2 ||/ ni2 || 2 < 2" +1 ||A: i ||™ +1 || Kt || 2 
As for the time-derivative of f n +i,i, we find 

f n+h i(x,z) =~\ J dy(h( x iy) (U,2(z,y) + fn,z(y,z)) + h(x,y) (j n ,i(z,y) + U,2(y,z)j 

+ h(y, z) (/„, 2 (x, y) + /„, 2 (y, x)) + fc 4 (y, z) (f n ,2(x, y) + /„, 2 (y, sen ) 
and similarly for f n +i,2- Hence, we find 

||/n+l J l||2<2f||^|| 2 ||/ nj2 || 2 + ||^|| 2 ||/ ni2 || 2 



< 2 (21^||^||^||I +4 n ||fc t ||^(||fc t || 2 ||fc t || 2 + 11^111 
n+l , /inMlj,. nniir, i|2 , _ ,n || u lin+1 



<(2^+4 ra )||fc t ||^||fe 4 ||^ + 2.4™||fc t ||r i H^ 
^^HWfllfellallMa + llMl 
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proving (|1 17[) for i = 1. The same bound for i = 2 and the second bound in (|118p for i = 1,2 
can be proven similarly. 

If n is odd, we have, using again the canonical commutation relations, 

ad« +1 (S) 

= ~ J dxdy (kt(x, y)a* x a* y - k t (x, y)a x a y ^j , i /" dxdy y)a*a y + /„, 2 (x, y)a x a* y ) 

= \ j dxdz(a* x a* z f n+1A (x,z) + a x a z f n+lj2 (x,z)) 
where 

/ n +i.i(z,z) = - / dyk t (x,y)(f nA (z,y) + f n ^(y,z)) , 

(121) 



/n+i,2(x,z) = - / dykt(x,y) (f n ,i(y,z) + f n ,2(z,y)) ■ 

The bounds (fTTHD . ([TT7jl . (fTTHjl follow as above. □ 

Using Lemma 16.91 we obtain the following properties of (dtT*)T. 
Proposition 6.10. There exists a constant C > with 

±(id t T*)T< Ce K ^(Af + l), 
± [Af, (id t T*)T] < Ce m (TV + 1), (122) 
±d t [(id t T*)T] <Ce m (Af + l). 

Proof. We observe first of all that, for all fx, f 2 G L 2 (R 3 x M 3 ) with J dx |/2(x,x)| < oo, we 
have 



ip, j dxdy (fi(x, y)a* x a* y + f 2 (x, y)a x a y ) i\> 

and 

{ip, / dxdy {fi(x, y)a* x a v + f 2 (x, y)a x a* y )i) 



<(||/i||2 + ||/2|| 2 )<^(Ar+l)^> (123) 



(124) 

<(||/l|| 2 +||/2|| 2 )<^^>+ / \f 2 {x,x)\dxU\\ 2 . 



In fact, (j!23|) follows because 

(V>, j dxdy (fi(x, y)a* x a* y + f 2 (x, y)a x a y )^ 

< j dx {\\aM\\o*{h{x,.))n + h*(f2(x, .))^|| \\a x 
<[|(JV + 1)^11 / efadl/i^Olla + ll/a^OllaJK 

<(||/l||2+||/2||2)||(AA+l) 1 /^|| 2 . 
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Eq. (|124p can be proven similarly. Combining the last estimates with Lemma 16,91 and with 
jnU), we find 



n>0 ^ 

^ E7^TTT! a l^ill2 + ll/n, 2 || 2 )||(AA + l 

">0 ^ 
n>l v 7 v,y 



n>0 

+ ^(2n)! 1/ da; l/2n-i,2(x,x)|(ix 

<cE»iifeibiK^ + i)'/w + i:«ii4ii 2 

<Ce 2 ll**ll a ||fc t || 2 ||(JS/'+l) 1/ Vl| 2 
^Ce^lKAA+l) 1 / 2 ^! 2 

using also Lemma l3,4i The second inequality in (|122p follows similarly because, essentially, 
the only consequence of taking the commutator with J\f is to eliminate the terms ad^(5) 
for all odd n. Also the third bound in (|122p can be proven analogously, taking the time 
derivative of the expressions for adg(.B) given in fj 1 1 5 j) . using the bounds for ||/ n ,j||2 in (|117p 
and (|118p and, finally, using the estimate for \\kt\\2 proven in Lemma 13, 41 □ 

6.6 Proof of Theorem 13.51 

In this section, we combine the results of the previous subsection, to obtain a proof of 
Theorem 13.51 From ([77|) and Propositions 16.11 16.31 16. 5( 16.61 it follows that 



C N (t) = C N (t) + K + ^ J dxdyN 3 V(N(x - y))a* x a* y a y a x 
"iV 3 / dxdy(Aw)(N(x-y))<pi N \x)<p { t N) (y)a* x a* y 



+ 

1 

2 

+ £(t) 



v 

+ - I dxdyN 3 V(N(x-y))(l-w(N(x-y)))^ N \x)4 N) (y)a x a y +h.i, 



(125) 



where the constant CV(i) is defined in (|58p and the error £ (£) is such that, for every 5 > 
there exists C5 > with 

±£{t) <5^IC + J dxdyN 2 V(N(x - y))a* x a* y a y a x ^j + C s ^- + C^ 1 * 1 (Af+1), 

± [Af,£(i)] < <5 (V + y dxdyiV 2 y(iV(x - y))a* x a* y a y a x ^ + C a ^r + C^ 1 * 1 (W + 1) , 

±£ (t) < <5 ^/C + y dxdyN 2 V(N(x - y))a* x a* y a y a x ^j + C^ 1 * 1 +AA+ 1^ . 

(126) 

In deriving (jl25f) . we made use of the crucial cancellation between the linear contributions in 
([77]) and the linear terms in ()102p . Next, we notice another crucial cancellation. The terms 
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on the third and fourth line in (|125p can be written as 
iV 3 J dxdy a* a* 

since / = 1 — w is a solution of the zero-energy scattering equation (—A + (1/2) V)/ = 0. 
We conclude that 

C N (t) = C N (t) + K + ^ j dxdyN 3 V(N(x - y))a* x a* y a y a x + 8{t) (127) 

where the error £(t) satisfies (|126p . Then (|60p follows from the first bound in (|126p . taking 
5 = 1/2. Also (|62p and (|63p follow from the second and third bounds in (|126p . since both fC 
and the quartic term on the r.h.s. of ()127p commute with N and are time-independent. 
This concludes the proof of Theorem 13.51 



+ 5 V|(1 



w 



(N(x-y))4 N) (x)4 N) (y) = 



A Properties of the solution of the Gross-Pitaevskii equation 

Proof of Proposition \3.1[ (i) This part of the proposition is standard. One proves first local 
well-posedness of the two equations in if 1 (lR 3 ). The time of existence depends only on the 
i? 1 -norm of the initial data. Since V, f > 0, the iJ 1 -norm is bounded by the energy, which 
is conserved. Hence one obtains global existence and a uniform bound on the i7 1 -norm. 

(ii) Also this part is rather standard, but since the non-linearity in (|50p depends on N, 
and we need bounds uniform in N, we sketch the proof of the bound (|5ip for ||</>t \\h u (the 
bound for ||</?t||# n can be proven analogously). We present the proof for the case t > 0. We 
claim, first of all, that there exists T > depending only on ||</?||#i and n G N such that 

sup H^^ll^n < 2||<po||ff* + sup H^f^llfcj-n-i. (128) 
te[o,T] te[o,T] 

Introducing the short-hand notation Un{x) = N 3 V(Nx)f(Nx), we write the solution tpf 1 ' 
of J5JU as 

Jo 

Differentiating this equation w.r.t. the spatial variables we find that 
9 * lp f) =e itA d a lp _ i f dse^-^Y^Y. (fl f)^*^^))^^. 

Here a is a three-dimensional multi-index of non-negative integers, with \a\ < n. 

The L£°([0, T], L^)-norm of the above expression can be controlled using Strichartz esti- 
mates for the free Schrodinger evolution e ztA (see Theorem 1.2]). We find 

\\o </>(.) \\L?Ll 

< ir*. + EE(«) (!) ii P» - ar^-^))if-^ 

< \\p<pw„+tw y: e (fi i ) -p ii * (d^d^f ^a-^x^. 
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By Holder and Young inequality, we find 



p<au</3 *e[o,T] 



for pi,P2,P3 > 1 with p x 1 + p 2 1 + £> 3 1 = 5/6. It is important to note that the indices 
(piiP2-,Pz) can be chosen differently for each term in the summation. In some of the terms 
with \a\ = n, all n derivatives hit the same (ft • Since 1/6 + 1/6 + 1/2 = 5/6, these terms 
can be bounded by 

ll^lllall^lU^CII^H^ (129) 

for C depending only on ||</?||#i (recall here that the if 1 -norm is bounded uniformly in t, by 
part (i)). In some of the other terms, one (p[ N ^ has n — 1 derivatives, one has at most one 
derivative and the last one has no derivatives. Since Hd 7 ^^ 1 1 < \\ip[ \\h u , if |t| < n — 1, 
these terms are bounded by the r.h.s. of (|129p . In all other terms, the three copies of <ft 
have at most n — 2 derivatives. These terms are bounded by 

for all multiindices 71, 72, 73 with | — | < n — 2, for i = 1,2, 3. We conclude that 

\\&*<Pn ) \\L°°L* < \\d a V \\ L 2+CT 1 ' 2 sup ||<^ iV) ||^ + CT 1 ' 2 sup ||¥>f°llfl»-i- 
u te[o,T] te[o,T] 

Summing over all a with \a\ < n, we find 

sup y { t N) \\ H n < M^+CT 1 / 2 sup \W[ N) \\ Hn +CT 1 ' 2 sup H^ll^.!. 
te[o,T] te[o,T] te[o,T] 

Choosing T > so small that CT 1 / 2 < 1/2, we find (fT28l) . 

To show (|5ip . we iterate now ()128p . We proceed by induction over n. For n = 1, the claim 
follows from part (i). Suppose now that ||Vt^||#(n-i) — C n -\ exp(if n _i|i|), for constants 
C n _i,if n _i depending on H^llf^n-i) and, respectively, on ||<£>||#i. Let T be as in (|128]h For 
an arbitrary t > 0, there exists an integer j £ N such that (J — 1)T < t < jT. Then 

II ( N )\\ ^ II (AO II 

se[(i-l)TjT] 

<2||^J2 1)t ||h» + 2 sup H^ll^-i 

S £[(j-l)T,jT] 



s 0\\,J N ) II 1 or3 c 3X n _ijT 

< 2 l^ (i _i) T llH" +2C n _ 1 e 



Similarly we have 



Iterating j-times, we obtain 

£=0 
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for some constant C n depending on ||v||i? n an d K n depending only on ||</?||^i. 

(iii) From the modified Gross-Pitaevski equation (f5Q|) . letting Un{x) = N s V(Nx)f(Nx), 
we find 

I, .(Af) n , || (iV) || , ( TT | (A0|2\ (N) 

\m h<\m IU 2 + [ u N*m I )<Pt 



<rf°lla» + 



ft Woo 



u N *\4 N) \ 2 
< \\fi N) \\ H > + cii^iun^llill^lloo 

<C\\4 N) \\% 2 <Ce K W 

for a constant C depending only on ||^||h 2 an d ||C^v[|ij and for K > depending only on 
ll^Hij-i. Here we used part (ii). Applying a gradient to (f50j) . we find 



(TV) 



Vft 



(TV) 



(130) 



Clearly, ||VAyj^ [|2 < \\if[ \\h 3 - The second term on the first line is bounded in norm by 
\(U N * \<pi N) \ 2 )V V f\ < \\(U N * |^| a )||oo||V^[| 2 

^ll^llill^llLllv^lb^cii^ii^. 

The terms on the second line of (|130|) can be bounded similarly. From part (ii), we conclude 
that ||V^^||2 < Cexp(K\t\). Analogously, we can also show that || V 2 ^^ || 2 < Ce K ^. We 
conclude that H^ - ^ ||#-2 < Ce K ^. Finally, ([50]) implies 



(TV) 



+ 



(u N * (iip 



t ft 



ipl + [U N * \if\ 



m 



Plugging in the r.h.s. of (|50p for i(p[ N \ we arrive at 



-ft ' = AVf } - A (Vtv * l^l 2 )^) + (U N * l^l 2 ) V 



+ (U N * #)(-A^))^ + 2 [tfc * (\fi N) \ 2 (UN * l^l 2 ) 
+ * (-A^VrVf) + * IvT I 2 ) (-A^)- 



ft 



(TV) 



Proceeding similarly as above, we find that \\^ { t N) \\ 2 < C\\<p[ N) \\ Hi < Cexp(K\t\ 
(iv) Using ([32]) and 0, we find 

dt\\(p t - ff^Wl = -21m (ip t , (U N * \^ t N) \ 2 - Svraol^tl 2 )^^ 
= -21m /(ft, (U N * \p t \ 2 ~ 87ra |^| 2 )^ Ar) ) 



(131) 



21m i(f U U N * (\ipl 



ftnif^) 



The second term on the r.h.s. can be written as 



Im ( Vt , (U N * (|^| 2 - M 2 )) rf*>) = Im (<p t , (U N * (l^l 2 - |^| 2 )) (f t - 4 N) )) ■ 



52 



Hence, by Holder's and triangle's inequality, 

< 



< 



u N *(\A N) \ 2 -\<Pt\ 2 ))(vt-vl N) ) 

UNH\4 N) \ 2 -\^\%yt-^ t N) h 

<MilNloo[|^-^l|2||^| a -N 2 | 2 

< ll^v||l|MU (iMloo + H^Hoo) \Wt - ^111 

<c(M^ + 11^11^) ll^-^ll!- 

As for the first term on the r.h.s. of (|13ip . we find (since J U^{y)dy = 87rao) 
(p t , (U N * \(p t \ 2 -8TTa \(p t \ 2 )(pl N ^ 

dxlp t {x)ipf\x) j dyU N (y) (\ip t (x - y)\ 2 - \ip t (x)\ 2 ) 
dxdyU N (y) \(pt{x)\\ipf\x)\ \\tp t (x - y)\ 2 - \yt{x)\ 2 \ . 



(132) 



< 

Writing Un(x) = N s U(Nx), with U(x) = V{x)f(x) and changing integration variables, we 
find 

<pt, {U N * \(pt\ 2 -8Tra \(pt\ 2 )(pi N)S ) 

< J ^[/(^l^^ll^^ljl^^-y/iV)! 2 -!^^)! 2 ). 



Using 



\\y t (x-y/N)\ 2 -\ip t {x)f 



ds—\tp t {x- S y/N)\< 
i ds 

<2|y|iV- 1 f ds\V<pt(x-8y/N)\\<pt(x-sy/N)\ 
Jo 



we conclude that 

[<pt,(U N * \(f t \ 2 - 87rao|<pt| 2 )<p^ 

< 2.V-' I dxdy f\sU(y)\y\\^ t {x)\\^\x)\\V Vt (x-sy/N)\\^ t {x-sy/N^ 

" x - (133) 



< 2iV- 1 ||^|| 2 J dxdy j dsU{y)\y\(\v[ N \x)\ 2 + \Vvt{x-sy/N)\' 



J tlll 



2 

H' 2 



where the constant C depends on j dyU(y)\y\ and on ||y||#i. Inserting (|133|) and (|132p into 
(|13ip . and using the estimate from part (ii) for ||<^||#2, we find 

ft||^-^||3<Ce*l*l||^-^||i + ^«l 

The claim now follows from GronwalFs inequality, since (ft=o = ^[=0- '— ' 
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B Properties of the kernel k t 

This section is devoted to the proof of Lemma 13.31 and Lemma [3? 



Proof of Lemma \3.3[ (i) We will make use of the bounds (|54p . The first bound implies 
immediately that 

^(^^l^mm^iVI^Cx)!!^^)!,^^!^^)!!^^)!) (134) 
and therefore, by Hardy's inequality 

||M£ <cj dxdy ] ^w{ N \ x )\^ ( ?\y)\ 2 < cil^ll^ll^lli < c. 

As for the gradient of kt, we have 

VM^V) = -N 2 Vw(N(x - y))4 N) (x)4 N) (y) - Nw(N(x - y))V^ 7V) (x)^ N) (y) 
and thus, from the second bound in ([54"]) . 

\\V 1 k t \\l<cj dxdy {N2{x + 1)2 \<pi N Hx)\ 2 \4 N \y)\ 2 

- CN I dx wwTTr yitN)rH1 - CN 

where we used Young and then Sobolev inequalities. Next we compute 
V 1 (k t k t )(x,y) =V X J dzk t (x,z)k t (z,y) 

<pi N \x)^ t N \y) I dzN 2 w(N(x-z)) W (N(z-y))\4 N \z)\ 2 



V^ t N \x)^ N \y) J dzN 2 w(N(x-z))w(N(z-y))\ V ( t N) (z)\ 2 
+ 4 N) (x)4 N) (y) [ dzN 3 Vw(N(x-z))w(N(z-y))\4 N) (z)\ 2 . 



Using ([54]) . we find 



\\Vl(k t k t )\\i <C dxdydz x dz 2 



+ C j dxdydz\dz 2 
< C / dxdydz\dz 2 



\x - z\\\z\ - y\\x - z 2 \\z 2 - y\ 

\^\x)\ 2 \^\y)\ 2 \vf\zi)\ 2 \vF\^ 



+ C I dxdydz\dz 2 



\x - z\\ 2 \zi -y\\x- z 2 \ 2 \z 2 - y\ 
\x - zi\ 2 \z 2 - y\ 2 

\^\ X )?\^\y)\ 2 W?\^)?W { r\^)? 



\x — z\\ 2 \z\ — y\ 2 \x — z 2 \ 

<c\\^\\Uw\ N) \\l 
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(ii) The pointwise bound for kt(x,y) follows directly from (|54|) . as noticed in (|134p . To 
bound \r(kt)(x,y)\, we observe that, by Holder's inequality, (f5"ij) and part (i), 

\(ktk t ) n kt(x,y)\ 

= \hk t {ktkt) n ~ l h{x,y)\ 
= \^\x)\\^\y)\ 

d Zl dz 2 N 2 w(N(x - z l )) W (N{z 2 -y))^f\z 1 )^\z 2 )k t {k t k t ) n - l {z 1 ,z 2 ) 

<\^ N \x)\W[ N \y)\\MhW- l h 

X ( I dz 1 N 2 w(N(x-z 1 )) 2 \i P ' t N \z 1 )\ 2 [ dz 2 N 2 w{N{z 2 -y)) 2 \^ 1 N \z 2 )\ 2 " 



Thus 

oo ^ 

\r(h)(x,y)\ < \(k t k t rk(x,y)\ < C\<p™ {x)\\ V ™ ^Wb. 

n=l 

The pointwise estimate for p(kt)(x,y) can be proven similarly. This completes the proof of 
part (ii). Part (hi) follows easily from the pointwise bounds in part (ii). □ 

Proof of Lemma \3.4\ In the following proof we will use the bounds ||^ \\h 2 i Wfit lb < 
(j e K\t\ f rom Proposition 13.11 
(i) From ([56]), we find 

INI 2 < 4 / dxd y] ^\^ N \x)\ 2 \^\y)\ 2 < CII^Hlll^ll^ < Ce^ (135) 

by Hardy's inequality, and from Proposition 13.11 part (hi). Similarly, 

ll^lb^Cll^lb + Cll^lbllV^lb^Ce^l 
by Proposition 13.11 Writing p(k t ) = ]Cn>o(^*^*) n /(2 n 0' we find immediately 

^ Z E (^yyHNblNir' 1 < \\hhe llhh < Ce K ^ 



n>0 



applying f|135j) . The bound for r(kt) can be proven analogously, 
(ii) From [231 part (v), we have 



[|Vip(fct)|| 2 < Ce^ [\\k t \\ 2 WVxihktjh + C||Vi(fctfc t )[| 2 + C[|Vi(fctfci)| 
We are left with the task of estimating || Vi(fct&t)||2 and ||Vi(fctfct)||2- We start by applying 
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the product rule: 

|2 



dxdy 



V x I dz (Nw(N(x - z))^ i t N \x)Tp < t N) (z)+Nw(N(x - z)) Tpf\x)Tpf\z) 
xNw(N(z-y))^ t N) (z)<pi N \y) 



< 4 / dxdy 
+ 4 J dxdy 
+ 4 / dxdy 



dzN 2 X7w(N{z - x))<? { t N \x)\<pf>(z)\ 2 ipf>(y)Nw(N(y - z)) 



dzNw(N(z - x))V^ N \x)\cp[ N \z)\ 2 ip[ N) (y)Nw(N( y - z)) 



dzN 2 Vw(N(z - x))Tp { t N] \x)i$*\z) 



(136) 
(137) 

(138) 
(139) 

Next, we estimate the four terms on the r.h.s. of the last equation. For the summands (|137|) 
and (|139p we use that, from (f54"|) . Nw(Nx) < C|x| _1 . Applying Hardy's inequality, both 
terms are bounded by C\\V(p { t N) \\l < Cexp(K\t\), using Proposition 13.11 Since, again by 
(JMD, N 2 Vw(Nx) < C\x\~ 2 , the contribution (fT36|) is bounded by 

2 



+ 4 / dxdy 



x4 N \ Z )4 N \ y )Nw(N(y-z)) 



dzNw(N{z - y))Vvl N \x)^ N \z)^ N) (z)cp < t N \y)Nw(N(y - z)) 



C I dxdy 



dz- 



1 



-\^\x)\\^\z)\ 2 \^\y)\ 



x — z\ 2 \z — y\ 

\^\-)\ 2 \^\v)?W\ N \^)?WT\^ 



C I dxdydz\dz 2 



\z\ - y\\z 2 -y\\x- zi| 2 |x - z 2 \ 2 



C I dx\v { t N) {x)\' z I dz x dz 



\^\z,)?\^\z^ 
\x — zi| 2 |x — z 2 \ 2 



dy- 



\z\ - y\\z 2 - y\ 

<C||^||i||^||^<Ce^l 

by Proposition 13.11 Analogously, we can also bound the contribution (|138p . This shows 
the bound for ||Vip(fct)||2. The bounds for ||V 2 p(^)||2, ||V2p(^)||2, ||V 2 p(^)||2 are proven 
similarly. 

(iii) From (|56p . using Nw(Nx) < C|x| _1 , we find immediately that 

supll^aOHa < C (llV^lbll^lU + ||V^|| 2 ||^||oo) < Ce K ® 
by Proposition 13.11 To show the bound for p(kt), we observe that 

p(k t )(x,y) = d t 77 ^-tt / dzidz 2 k t (x,z 1 )(k t kt)( n ~ 1 \zi,z 2 )k t (z 2 ,y) 



E 

n>0 



(2n)! 



dz x dz 2 k t (x, zi)(/c t /c t ) (n 1 \z 1 ,z 2 )k t (z 2 ,y) 

+ J dz 1 dz 2 k t {x,z 1 )(d t (kth) {:n '~ V) ){zi, z 2 )k t {z 2 , y) 
+ / dzidz 2 k t (x,zi)(k t k t ) ( - n ~ 1) (z 1 ,z 2 )k t (z 2 ,y) . 



(140) 
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The first term in the parenthesis can be bounded in absolute value by 

dz 1 dz 2 h(x, z 1 )(k t k t ) in ~ 1 Xz 1 , z 2 )k t (z 2 ,y) 
< C [ d Zl dz 2 r-J— T (^(xJH^Czi)! + irf^ll^f } (^)| 



X — Z\\ 

x \(k t k t )^- 1 \zi,z 2 )\ -^—A^ (z^W^iy)] 

\y - z 2\ 

<C\4 N \x)\\<pl N \y)\ I dz x dz 2 _ \ _ \^\z l )\\^\z 2 )\\{kM n - l \z ll z 2 )\ 
J i'-e z\\\y z 2 \ 

+ C\^ t N \x)\\4 N \y)\ I dz x dz 2 _ \ _ - \^f\z l )\\^\z 2 )\\{k t ht- 1 \zuz 2 )\ 

j i-c ^iiiy z 2 1 
<c\\Chhr-% (iv-^^-)!!^^^^)!!!^^^^!!!^^^^! -k I v^^^ C^) 1 1 C-/) 1 1 1 ^v-^^ 1 1 2 ) - 

The last term in the parenthesis on the r.h.s. of (|14(jp can be bounded analogously. The 
middle term, on the other hand is bounded in absolute value by 

dz x dz 2 k t (x, z 1 )(d t (k t k t )( n ~ 1 ' ) )(zi,z 2 )k t (z 2 ,y) 

<C\^ N \x)\\^ N \y)\ [ d^d^ 1 ^^)''^^^ 1 ^^^)"- 1 ^,^)! 
J \x-zi\\y-z 2 \ 

< C\\d t (k t k t r-% \\V4 N) \\ 2 2 \^ t N \x)\\4 N) (y)\ 

Inserting the last bounds in ()140p . we find 

\p(k t )(x,y)\ < Ce K \%^ (U N \x)\ + \<pl N \x)\)(\^ N \ y )\ + |# } (y)|). 
Integrating over x and taking the supremum over y gives, as before using ()5ip . 

S up\\p(k t )(.,y)\\ 2 <Ce K W. 

y 

The bound for r(kt) can be proven analogously. Combining the bound for r(kt) with the one 
for kt, we also obtain the bound for sh(fc(). □ 

C Convergence for TV-particle states 



In this section, we show how the result of Theorem 11.11 stated there for initial data of the 
form W(\^Nip)T(ko)ip can be extended to a certain class of data with number of particles 
fixed to N. 

Theorem C.l. Let (p G H 4 (R 3 ) and suppose ip 6 J- with \\ip\\j^ = 1 is such that 

^(^+m + Hn)A <c (mi) 

for a constant C > 0. Let Pjy denote the projection onto the N -particle sector of the Fock 
space and assume that 

\\P N W{y/Nip)T(k )ip\\ > CN- 1 ^ (142) 
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for all N 6 N large enough. We consider the time evolution 

,-m N t PNW(VN<p)T(k )iP 



P N W{VNv)T{k 



and we denote by the one-particle reduced density associated with the N -particle state 
ipNj ■ Then there exist constants C, c\ , C2 > with 



Tr 



C exp(ciexp(c 2 |t|)) 

ATl/4 



for all t € K and all N large enough. Here (ft denotes the solution of the time- dependent 
Gross-Pitaevskii equation with initial data (ft=o = <*P- 

Remarks. 

(i) If we relax (fT32"|) to the weaker condition \\P N W(VN(p)T(k )^\\ > CN~ a , for some 
1/4 < a < 1/2, the proof below still implies the convergence 7^ — > \(pt)((ft\ but only 
with the slower rate N~ l l 2+a . 

(ii) The assumption (|142p and its weaker versions mentioned in the previous remark are 
very reasonable; let us explain why. The expected number of particles in the Fock state 
W(y/N(p)T(k )i/) is given by 



W(VN(p)T(ko)ip,MW(VNtp)T(ko)i/>) ^ 
= N + VN{T(k )^, H<p)T(k )iP) + (T(k )^,MT(k ^) 

with the notation <ft((f) = a(ip) + a*((p). From Lemma [2. II Lemma [4. 3 1 and the assump- 
tion (|14ip we conclude that there exists a constant C > with 



N - CN 1/2 < (w(VN<p)T(ko)^,MW(VN(fi)T(k )^J <N + CN 1/2 . 
The expectation of Af 2 , on the other hand, is given by 
(W{VN^)T{k )^,Af 2 W(VNp)T(k )i(j) 



= (T(k ^, (Af + VN<f>(ip) + N) 2 T(k )^) 

= N 2 + 2N 3 / 2 (T(k )^, <t>((p)T(ko)iP) + 2N(T(k ^,MT(k )iP) 

+ N(T(k )i>, <%) 2 r( W) + VN(T(k )4>, (Af^(ip) + 4>{<p)N) r(fco)V) 

+ (T(k )4>,Af 2 T(k )t/j). 

Subtracting the square of ()143p , and applying again Lemma I2.1| Lemma 14.31 and the 
assumption (|14ip . we estimate the variance of the number of particles in the state 
W(VNip)T(ko)*l> by 



W{VN(p)T(ko)ip, (Af - {N)yW(VNtp)T(ko)il>) < CN 



for an appropriate constant C > 0. Here (Af) is a shorthand notation for the l.h.s. of 
dUSl). We conclude that 

W(VNip)T(k )ip, 1 (\Af -{Af}\> KVN^ W(VNip)T(k )ip) < CK~ 2 . 
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Choosing K > sufficiently large, we find 

(w(VNip)T(k )*p, 1 (\Af - (A/"> I < ^Viv) W(y/N<p)T(ka)if>} > 1/2. 
From (|143p . adjusting the value of K, we obtain 

V(v^^)r(fc )^, 1 (|AA - AT | < W(y/N<p)T(ko)il>) > 1/2. 



2 



This means that 

n+kVn 

\Pj W(VN(p)T{k )ip " > 1/2. 



The average value of ||Pj WfyASO^Oo)^!! for 3 between N - K^N and N + iWiV 
is therefore larger or equal to iV -1 / 2 , in accordance with the assumption (I142p . In fact, 
this argument shows that for every N there exists an M G [iV — K^fN,N + iT\/iV] 
with ||P M W(V3Vy>)T(fco)V>|| > iV- 1 / 4 . Letting 



WN,M,t = e 



|P M W(VNip)T(k 

j 

can show, similarly to Theorem lC.il that 



and denoting by ^^ Mt the one-particle reduced density associated with ViV,M,tj one 



Tr 



C 1 exp(ci exp(c 2 |t|)) 

ATl/4 



The fact that the number of particles M does not exactly match the parameter N 
entering the Hamiltonian and the Weyl operator W(yN<p) does not affect the analysis 
in any substantial way, since \M-N\ <CN l l 2 <CiV. 

Proof of Theorem \C.1\ We write the integral kernel of as 



TnA x >V) 



1 



N\\P N W(VN<p)T(k 

x (e- iHNt P N W(VN(p)T(ko)ip, a* y a x e- mNt P N W(VN<p)T(k )^ 



1 



JV||PjvW(vJVfyOT(Jfet) 

x (e- mNt P N W(VNv)T{ko)ip, W(VN^ { t N) ) (a* y + v^V#°(y^ 

x (a x + \/iV^ (Ar) (x)) W(\/iV^ 7V) )e-^ JV *Ty(v / iV^)T(/t )v) • 
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Hence, we find 

ft{x) 



+ 



N\\P N W(VNtp)T(k )i>\\2 

x (e- mNt P N W{^N^)T{k )^, W{^^ ) )alW{^^ ) )e^ mNt W{VNip)T(k G )^ 



'N\\P N W{^N^)T{k Q ) 
x (e- mNt P N W{^N^)T{k )^, W{^^ ] )a x W{^ip\ N) )e- mNt W{VNip)T{k^ 

xa*a x W{VNvl N) )e- mNt W(VNv)T{k )^ 



Therefore, for any compact operator J on L (Ur) we find 

Tr J (tS - \Vt){Vt\) 

- inNt P N W{VNcp)T(k )i(j, W(VNip { t N) 



N\\P N W(VN^)T(k )^\\ 2 

x (a(Jip t ) +a*(Ji Pt ))W(VNip ( t N) e- mNt W(VNtp)T(k )iP 



xdT(J) W(^/Nip { t N) e- mNt W(VNip)T(ko)^ 

where dT(J) denotes the second quantization of the operator J, defined by (dT(J)t/j)^ = 
Ya=i J®i>^ for every ip = {ip^} n >o G J 7 (here J® denotes the operator acting as J on the 
z-th particle and as the identity on the other (n— 1) particles). Since ||<ir( J)ip\\ < ||J| 
we find, applying Lemma |2.1| 



Tr J (7$ - i^x<pti) 



where ||J|| denotes the operator norm of J. From Lemma 14.31 recalling the definition ([5 
of the fluctuation dynamics, we find 



T*j( f r$t-\<pt)(<f>t\) 



< —= ill \\M 1/2 T*{k t )W*{VN^ ( t N) )e-' mNt W{VN^)T{k Q 



N\\P N W(VNp)T(k 



+ 



1 



N\\P N W{VNip)T(k 



\\Af T* (kt) W*(VN(pl N) e- mNt W(VN(p)T(k 



N\\P N W(VN<p)T(k )i(j\\ N\\P N W(VN(p)T(k 

GO 



Using Proposition 14.21 we conclude that 



J [TnI ~ \<Pt)(<Pt 



< 



C\\J\\\\{M+lfl 2 U{t; 



N\\P N W(VNif)T(k 
C\\(M+m 



+ 



CWiN+l) 1 ' 2 ' 



N\\P N W(VNcp)T(k 



N\\P N W(VN(p)T{k ' 
From the assumptions ()14ip and ()142p . we obtain 

7$ - \<pt){<Pt\)\ < SllCAf + i) 1/2 ^;0)VII + 



Tr J 

Finally, Theorem 14.11 implies that 



C 



ATI/4" 



Tr J {in\ - I I ) 



< 



C|| J\\ exp (ci exp(c2|t|)) 



(144) 



Since the Banach space C 1 (L 2 (M. 3 )) is the dual space to the space of compact operators, 
equipped with the operator norm, (|144p implies the claim. □ 
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